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Stratospheric Temperature Changes: Observations and Model 

Simulations 


Summary 

This paper reviews observations of stratospheric temperatures that have 
been made over a period of several decades. Those observed temperatures have 
been used to assess variations and trends in stratospheric temperatures. A wide 
range of observation datasets have been used, comprising measurements by 
radiosonde (1940s to the present), satellite (1979 - present), lidar (1979 - present) 
and rocketsonde (periods varying with location, but most terminating by about 
the mid-1990s). In addition, trends have also been assessed from meteorological 
analyses, based on radiosonde and/or satellite data, and products based on 
assimilating observations into a general circulation model. 

Radiosonde and satellite data indicate a cooling trend of the annual-mean 
lower stratosphere since about 1980. Over the period 1979-1994, the trend is 
0.6K/ decade. For the period prior to 1980, the radiosonde data exhibit a 
substantially weaker long-term cooling trend. In the northern hemisphere, the 
cooling trend is about 0.75K/ decade in the lower stratosphere, with a reduction 
in the cooling in mid-stratosphere (near 35 km), and increased cooling in the 
upper stratosphere (approximately 2 K per decade at 50 km). 

Model simulations indicate that the depletion of lower stratospheric ozone 
is the dominant factor in the observed lower stratospheric cooling. In the middle 
and upper stratosphere both the well-mixed greenhouse gases (such as CO 2 ) and 
ozone changes contribute in an important manner to the cooling. 



1. INTRODUCTION 


For at least a decade now, the investigation of trends in stratospheric temperatures has been 
recognized as an integral component of ozone change investigations (e.g., WMO reports since 
early 1980s). A comprehensive international scientific assessment of stratospheric temperature 
changes was undertaken in WMO (1988). Analyses of the then available datasets (rocketsonde, 
radiosonde, and satellite records) over the period 1979/80 to 1985/86 indicated that the ob- 
served temperature trend was inconsistent with the then apparent ozone losses inferred from 
SBUV data, but consistent with SAGE ozone changes. The largest cooling in the observed 
datasets was in the upper stratosphere, while the lower stratosphere had experienced no signif- 
icant cooling except in the tropics and Antarctica. It is interesting to note that the period ana- 
lyzed by the 1987 Temperature Trends Panel (1979-1986) was one when severe ozone losses 
were just beginning to be recognized in the Antarctic springtime lower stratosphere. It was also 
a period from sunspot maximum to sunspot minimum. 

Since the time of the 1988 WMO assessment, there has been an ever-growing impetus for 
observational and model investigations of the stratospheric temperature trends (WMO: 1990, 

1 992, 1 995). This has occurred owing to the secular increases in greenhouse gases and the now 
well-documented global and seasonal losses of stratospheric ozone, both of which have a sub- 
stantial impact on the stratospheric radiative-dynamical equilibrium. The availability of various 
temperature observations and the ever-increasing length of the data record have also been en- 
couraging factors. In addition, models have progressively acquired the capability to perform 
more realistic simulations of the stratosphere. This has provided a motivation for comparing 
model results with observations, and thereby search for causal explanation/s of the observed 
trends. The developments seen in modeling underscore the significance of the interactions be- 
tween radiation, dynamics and chemistry in the interpretation of linkages between changes in 
trace species and temperature trends. Temperature changes are also instrumental in the micro- 
physical -chemical processes of importance in the stratosphere (see Chapter 7). 

The assessment of stratospheric temperature trends is now regarded as a high priority in climate 
change research inasmuch as it has been shown to be a key entity in the detection and attribution 
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of the observed vertical profile of temperature changes in the Earth’s atmosphere (Hansen et al., 
1995; Santer et al., 1996; Tett et al., 1996). Indeed, the subject of trends in the stratospheric tem- 
peratures is of crucial importance to the Intergovernmental Panel on Climate Change assessment 
(IPCC, 1996), and constitutes a significant scientific input into policy decisions. 

An excellent perspective of the evolution in the state-of-the-science of the stratospheric 
temperature trends can be found in the WMO Ozone assessment reports - WMO (1986, 1988, 
1990, 1994). We summarize here the principal results concerning stratospheric temperature 
trends from these WMO reports. In general, the successive assessments since the WMO (1986) 
and WMO (1988) reports have traced the evolution of the state-of-the-science on both the ob- 
servation and model simulation fronts. On the observational side, WMO has reported on avail- 
able temperature trends from various kinds of instruments - radiosonde, rocketsonde, satellite 
and lidar. On the modeling side, since the 1986 report, WMO has reported on model investiga- 
tions that illustrate the role of greenhouse gases and aerosols in the thermal structure of the 
stratosphere, and the effects due to changes in their concentrations upon stratospheric temper- 
ature trends. WMO ( 1 990) began to recognize, from observational and modeling standpoints, 
the substantial lower stratospheric cooling occurring during springtime in the Antarctic as a 
consequence of the large ozone depletion. Low and middle latitude lower stratosphere were in- 
ferred to have a cooling of less than 0.4K/ decade over the prior 20 years. The upper 
stratosphere was estimated to have cooled by 1.5 +/- IK between 1985/86 and 1979/80. 

WMO (1992) reported that, based on radiosonde analyses, a global-mean lower strato- 
spheric cooling of -0.3K had occurred over the previous 2-3 decades. Model calculations indi- 
cated that the observed ozone losses had the potential to yield substantial cooling of the global 
lower stratosphere. At 44N during summer, a cooling of the upper stratosphere (-1K/ decade 
at -35 km.) and mesosphere (-4K/ decade) was also reported. However, in general, for assess- 
ment purposes, the global stratospheric temperature record and understanding of temperature 
changes were found to be not as sound as those related to ozone changes. 

The WMO (1995) assessment discussed the observational and modeling efforts through 
1994, focussing entirely on trends in lower stratospheric temperatures. This assessment con- 
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eluded, based on radiosonde and satellite microwave observations, that there were short-term 
variations superposed on the long-term trends. A contributing factor to the former was strato- 
spheric aerosol increases following the El Chichon and Pinatubo volcanic eruptions which re- 
sulted in an increase of the global stratospheric temperature. These transient warmings posed 
a complication when analyzing the long-term trends and deducing their attribution. The long- 
term trends from the radiosonde and MSU satellite data indicated a cooling of 0.25 to 0.4 YU 
decade since the late 1970s, with suggestions of an acceleration of the cooling during the 1980s. 
The global cooling of the lower stratosphere suggested by the observations was reproduced rea- 
sonably well by models considering the observed decreases of ozone in the lower stratosphere. 
For altitudes above the lower stratosphere, a clear conclusion concerning trends could not be 
made. The 1992 and 1995 assessments laid the basis for the conclusion that the observed trends 
m the lower stratosphere during the 1980s were largely attributable to halocarbon-induced 
ozone losses. 

In this review, we extend the evaluations of the aforementioned WMO assessments by fo- 
cussing on the decadal-scale lower-to-upper stratospheric temperature trends arising out of ob- 
servational (section 2) and model simulations (section 3) analyses. Towards this end, we utilize 
the results available in the recently published literature. The temperature observations consid- 
ered span at least 10 years, the period considered for evaluation is typically at least 15 years. 
The trend estimates discussed here include (i) a long-term period that spans two decades or 
more, (ii) the period since 1979 and extending to either 1994 or to the present (i.e., up to 1998), 
and (iii) the period extending from 1979 to about early 1990s. The last-mentioned period is that 
for which several model simulations have been compared with observations. In section 4, we 
use the results from sections 2 and 3 to investigate the extent to which the observed temperature 
trends can be attributed to changes in the concentrations of radiatively-active species. The dis- 
cussions and findings here also constitute Chapter 5 of the WMO (1998) Ozone assessment re- 
port. 


2. OBSERVATIONS 
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The types of observational data available for investigation into stratospheric temperature trends 
are diverse. They differ in type of measurement, length of time period and space-time sampling. 
There have been several investigations of trends that have considered varying time spans with the 
different available datasets, as will be discussed shortly. In recent years, the World Climate Re- 
search Program’s SPARC-STTA (Stratospheric Processes and their Role in Climate - Stratospheric 
Temperature Trends Assessment) group has initiated a project to bring together various datasets 
covering the period 1966-1994, and to intercompare the resulting global stratospheric temperature 
trends. SPARC- STTA chose two different time periods to examine the trends based on the avail- 
ability of the data viz., 1979-1994 and -1965-1994. The shorter period coincides with the period 
when severe global ozone losses have been detected and also coincides with the period of global 
satellite observations. The second period is a longer one for which radiosonde (and a few rocket) 
datasets are available. 

The updated datasets made available to and employed by the SPARC-STTA for the analyses 
are shown (with the exception of the rocketsondes) in Table 5-1 along with their respective lati- 
tudes, altitudes and periods of coverage. Additionally, independent of the STTA activity, some in- 
vestigations [Dunkerton et al. (1998), Keckhut et al. (1998), Komuro (1989), Golitsyn et al. (1996) 
and Lysenko et al. (1997)] have analyzed trends from rocketsonde observations made at a few geo- 
graphical locations and over specific time periods (see Table 5-2). We utilize these datasets in some 
of the presentations to follow. It is convenient to group the currently known datasets as follows: 

Ground-based instruments : radiosonde, rocketsonde, and lidar; 

Satellite instruments : microwave and infrared sounders; and 

Analyses', employing data from one or both of the above instrument types, without/ with a nu- 
merical model. 

The datasets indicated in Table (5-j)are a collection of monthly-mean, zonal-mean temperature 
time series. All but one of these datasets cover the years 1979-1994, and some extend further back 



in time. The pressure-altitude levels of the datasets vary, but overall they cover the range 100 to 0.4 
hPa (approximately 16-55 km.). Most datasets provide temperatures at specific pressure levels, but 
some provide data as mean temperatures representative of various pressure-layers. The instrumen- 
tal records from radiosondes, rocketsondes, lidar and satellite (MSU and SSU) are virtually inde- 
pendent of each other. General characteristics of the different datasets are discussed below (see also 
WMO: 1988, 1990). 

Radiosonde 

Radiosonde data are available dating back to approximately the early 1940s. Although the 
sonde data do not cover the entire globe, there have been several well-documented efforts to use 
varied techniques m order to obtain the temperatures over the entire northern hemisphere or the 
global domains. The sonde data cover primarily the lower stratospheric region (approximately, 
pressures greater than 10 hPa). The geographical coverage is quite reasonable in the Northern 
Hemisphere (particularly midlatitudes) but is poor in the extremely high latitudes and tropics, and 
is seriously deficient in the Southern Hemisphere (Oort and Liu, 1993). 

As was the case at the time of the 1 988 assessment, two organizations monitor trends and 
variatiqns in lower stratospheric temperatures using radiosonde data alone. The “Berlin” group 
(e.g., Labitzke and van Loon, 1 995) prepares daily hand drawn stratospheric maps based on synoptic 
analyses of radiosonde data at 100, 50, 30, and, in some months, 10 hPa, beginning from 1964. The 
“Berlin” monthly dataset examined by SPARC-STTA is derived from these daily analyses. The 
NOAA Air Resources Laboratory (e.g , Angell, 1988) uses daily radiosonde soundings to calculate 
seasonal layer-mean “virtual temperature” anomalies from long-term means, and uses these to 
determine trends since 1958 in the 850-300, 300-100, and 100-50 hPa layers (lower stratosphere). 
[“Virtual temperature” is the temperature of dry air having the same pressure and density as the 
actual moist air. Virtual temperature always exceeds temperature, but the difference is negligible 
in the stratosphere (Elliott et al., 1994)]. The layer-mean virtual temperatures are determined from 
the geopotential heights of the layer endpoints. The Berlin analyses are of the Northern Hemisphere 
stratosphere and troposphere and are based on all available radiosonde data, whereas the Air 
Resources Laboratory monitors trends at 63 stations, in seven zonal bands covering the globe. 
Additionally, Angell (1991b) also monitors stratospheric temperature, particularly its response to 
volcanic eruptions, at four levels between 20 km. (50 hPa) and 3 1 km. (10 hPa) using a network of 
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12 stations ranging from 8S to 55N. The “Angell” data used here represents a subset. 

Extensive analyses of radiosonde temperature data by the NOAA Geophysical Fluid Dynamics 
Laboratory (GFDL) (Oort and Liu, 1993; labelled “Oort”) and the U.K. Meteorological Office’s 
(UKMO) Hadley Centre for Climate Prediction and Research (Parker et al., 1997) have also been 
used for quantifying stratospheric trends. The GFDL database consists of gridded global objective 
analyses based on monthly means derived from daily soundings for the period 1958-1989 for the 
tropospheric levels and the 100, 70, 50, and 30 hPa levels in the lower stratosphere. Layer-mean 
trends for the 100-50 hPa layer are based on temperature data at the 100, 70, and 50 hPa levels. A 
subset of this dataset is used here. The UKMO gridded dataset (“RAOB” or “UKRAOB”) for 
1958-1996 is based on monthly mean (CLIMAT TEMP) station reports, adjusted (using MSU 
Channel 4 data, discussed below, as a reference) to remove some time-varying biases since 1979 
for stations in Australia and New Zealand, and interpolated in some data- void regions. The “Russia” 
set consists of data from the high northern latitudes (70 and 80N; Koshelkov and Zakharov, 1998). 
Thus, the “Angell, Berlin, Oort, Russia and Raob” datasets used by SPARC-STTA are compilations 
of data from various radiosonde stations, grouped, interpolated and/or averaged in various ways to 
obtain monthly-mean and latitude-mean, pressure-level or vertical-average temperatures. 

Rocket and Lidar 

Rocket and lidar data cover the altitude range from about the middle into the upper stratosphere 
and mesosphere. Rocketsonde data are available through the early 1990s from some locations but 
the activity appears to be virtually terminated except in Japan (see Table 5-2). The lidar measure- 
ment, just like the rocketsondes, has a fine vertical resolution. Lidar measurements of stratospheric 
temperatures are available since 1979 from the Haute Provence Observatory (OHP) in southern 
France (44N, 6E). Specifically, the “lidar” (Table 5-1) temperatures observed at altitudes of 30 to 
90 km. are obtained from two lidar stations, with data interpolated to pressure levels (Keckhut et 
al., 1995). Several other lidar sites have initiated operations and could potentially contribute in fu- 
ture temperature trends assessments. 

MSU and SSU Satellites 

Satellite instruments of interest have become available since -1979 (Table 5-1). These fall into 
two categories - remotely sensing in the microwave (Spencer and Christy, 1993) and thermal infra- 
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red (Nash and Forrester, 1 986) wavelengths. The “MSU” Channel 4 dataset derives from the lower 
stratosphere channel (-150-50 hPa) of the Microwave Sounding Unit on NOAA polar operational 
satellites (Figure 1 , left panel). The “Nash” dataset consists of brightness temperatures from ob- 
served (25, 26 and 27) and derived (47X,_36X, 35X, 26X and 15X) channels of the Stratospheric 
Sounding Unit (SSU) and HIRS-2 instruments on these same satellites (Figure 1). One complica- 
tion with satellite data is the fact that there are discontinuities in the time series owing to the mea- 
surements being made by different satellites monitoring the stratosphere since 1979. Adjustments 
have been made in the “Nash” channel data provided to compensate for radiometric differences, 
tidal differences between spacecraft, long-term drift in the local time of measurements, and also 
for spectroscopic dnft in channels 26 and 27. Adjustments have also been made to MSU data (e.g., 
Christy et al., 1995). 

An important attribute of the satellite instruments is their global coverage. However, in contrast 
to the ground-based measurements e.g., the radiosonde, which performs measurements at specific 
pressure levels, the available satellite sensors have response functions that sense the signal from a 
wide range in altitude. The nadir satellite instruments ‘sense’ the emission originating from a layer 
of the atmosphere approximately 10-15 km thick. Figure 1 illustrates the weighting function for 
the MSU and SSU channels analyzed here, exhibiting the thick-layer nature of the measurements. 
For example, the emission for the microwave MSU Channel 4 comes from the -12-22 km. layer, 
while, for the thermal infrared SSU Channel 15X, it is from the -12-28 km. layer. In Table 5-1, a 
nominal center pressure of each satellite channel has been designated, but it is emphasized that the 
preponderance of energy comes from a -8-12 km. thick vertical layer centered around the con- 
cerned pressure level, as indicated in Figure 1 . A perspective into the global-mean anomalies of 
temperature at various stratospheric altitudes between 1979 and 1995 (deviations with respect to 
the mean over this period), as derived from the different SSU channels, can be obtained from Fig- 
ure 2. The MSU record has been discussed in WMO (1995). 

Analyzed datasets 

There are a number of datasets that involve some kind of analyses of the observations. They 
employ one or more types of observed data, together with the use of some mathematical technique 
and / or a general circulation assimilation model, to construct the global time series of the temper- 
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atures. They are, in essence, more of a derived dataset than the satellite- or the ground-based ones. 

The “CPC” (Climate Prediction Center, formerly Climate Analyses Center or “CAC”) and “UK- 
MO/SSUANAL” stratospheric analyses (Table 5-1) do not involve any numerical atmospheric cir- 

A * 

culation model. The “CPC” northern hemisphere 70, 50, 30 and 10 hPa analyses use radiosonde A * 

data. Both the “CPC” and “UKMO/SSUANAL” analyses (see also Swinbank and O’N eill, 1994) 1 

use TOVS (TIROS Operational Vertical Soundings) temperatures, which incorporate data from the 
SSU, HIRS-2, and MSU on the NOAA polar orbiting satellites. Adjustments based on rocketsonde 
(Finger et al., 1993) have been applied to the “CPC” 5, 2 and 1 hPa temperatures for the SPARC 
dataset. [However, both CPC data above 10 hPa altitude and UKMO/SSUANAL datasets may be j’ 
limited in scope for trend studies as they have not yet been adjusted for some of the known prob- j 
lems associated with SSU satellite retrievals from many different satellites (see section 2.5.2); j 

U ir 

hence, they are not considered in the intercomparisons reported here]. -U- . 

JU-i- x ^ ^ 
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The “Reanal” (viz., the US National Centers for Environmental Prediction or NCEP reanalyzed) 
and the “GSFC” (NASA) datasets are derived using numerical atmospheric general circulation 
models as part of the respective data assimilation systems. These analyses projects provide synoptic „ 
meteorological data extending over many years using an unchanged assimilation system. In general, - 1 ty ) 

analyzed datasets are dependent on the quality of the data sources such that a spurious trend in a 
data source could be inadvertently incorporated in the assimilation. Also, analyses do not 
necessarily account for longer-term calibration-related problems in the data. Further, the analyzed 
datasets may not contain adjustments for satellite data discontinuities (Santer et al., 1998). 

2.2 Summary of various radiosonde trends investigations 

Radiosonde data show a cooling of the lower stratosphere over the past several decades. Table 
5-3 summarizes published temperature trend estimates by various investigators, including those 
mentioned in 5.2.1. The data periods and the analyses techniques vary, as do the levels and layers 
analyzed. No attempt is made here to critically evaluate these diverse estimation techniques. The 
reported trends have all been converted to units of K per decade. Overall, the trends, based on areal 
averages and all seasons, are negative and range from 0 to several tenths of a degree per decade. 

The few studies with global coverage show more cooling of the Southern Hemisphere lower 
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stratosphere than the Northern Hemisphere. Large trends evaluated for the decade of the 1980s 
emphasize the period of ozone loss. Positive trends have been found at a few individual stations in 
the tropics by Reid et al. (1989) for the period 1966-1982, possibly due to the influence of the El 
Chichon volcano effects. [It may be noted that Labitzke and van Loon (1995) find positive trends 
(not listed in Table 5-3) at high and low latitudes for the month of January]. 

The sensitivity of trend estimates to the period of record considered is evident from the time 
series of global or hemispheric mean lower stratospheric temperature anomalies (Angell, 1988; 
Oort and Liu, 1993; Parker etal., 1997). These data (Figure 3(top); Angell, 1988, updated; Halpert 
and Bell, 1997) show relatively high temperatures (particularly in the southern hemisphere) during 
the early 1960’s, fairly steady temperatures till about 1981, and relatively low temperatures since 
about 1984, with episodic warmings associated with prominent volcanic eruptions. Figure 
3(bottom) shows the global temperature anomalies from the MSU satellite. The evolution of the 
anomalies is qualitatively similar to the radiosonde (Christy, 1995), including the warming in the 
wake of the El Chichon and Pinatubo eruptions (WMO, 1995), followed by a cooling to somewhat 
below the pre-eruption levels. The long-term cooling tendency of the global stratosphere is 
discernible in both datasets although the satellite data exhibits less interhemispheric difference. 

The Berlin analysis (Figure 4) shows that the radiosonde temperature time series for the 30 hPa 
region at the northern pole in July acquires a distinct downward trend when the 1955-1997 period 
is considered in contrast to the behavior for the 1955-1979 period. The trend estimates are seen to 
depend on the end years chosen. The summertime temperature decreases in the high northern 
latitudes have been more substantial and significant when the decade of 1980s and after is 
considered; note that this does not necessarily imply a sharp downward trend for the other months. 

A few studies have examined radiosonde observations of extreme temperatures in the lower 
stratosphere. At Sodankyla Finland, Taalas and Kyro (1994) found an increase in the frequency of 
occurrence of temperatures below 195K at 50 hPa during 1965-1992. At both 50 and 30 hPa, over 
the northern hemisphere (10-90N), Pawson and Naujokat (1997) found a decrease m the minimum 
and an increase in the maximum daily wintertime temperatures during 1965-1996. They also found 
an increase in the area with temperatures less than 195K, and suggested that extremely low 
temperatures appear to have occurred more frequently over the past 15 years. 
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2.3 Zonal, annual-mean trends 


2.3.1 TRENDS DETERMINATION 

There is a wide range in the numerical methods used in the literature to derive trends and 
their significance. Most studies are based on linear regression analyses, although details of the 
mathematical models and particularly aspects of the standard error estimates are different. Dif- 
ferences in details of the models include the method of fitting seasonal variability, the number 
and types of dynamical proxies included, and the method used to account for serial autocorre- 
lation of meteorological data (e.g., multiple linear regression analysis (MLRA) model of Keck- 
hut et al., 1995; model used by Randel and Cobb, 1994). 

The SPARC-STTA group calculated the temperature trends (K / decade) from each of the 
datasets using autoregressive time-series analyses (maximum likelihood estimation method; 
e.g., Efron, 1982). The methodology consists of fitting the time series of monthly mean values 
at each latitude with a constant and six variables (annual sine, annual cosine, semi-annual sine, 
semi-annual cosine, solar cycle and linear trend). The derived trend and standard error are the 
products of this computation. The t-test for significance at the 95% confidence level is met if 
the absolute value of the trend divided by the standard error estimate exceeds 2. The results 
from the statistical technique used by SPARC-STTA have been intercompared with other meth- 
ods employed in the literature (A. J Miller, personal communication) and found to yield similar 
trend estimates. It is cautioned, however, that the estimates of the statistical uncertainties could 
be more sensitive to details of the method than the trend results themselves, especially if the 
time series has lots of missing data. 

An important caveat to the interpretation of the significance of the datasets is that the time 
series analyzed below, in some instances, is only 15 years long or, in the case of the lengthier 
rocketsonde and radiosonde record, up to -30 years long. In this context, it must be noted that 
the low frequency variability in the stratosphere, especially at specific locations, is yet to be ful- 
ly ascertained and, as such, could have a bearing on the robustness of the derived trend values. 
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2.3.2 50 AND 100 hPa 


Figure 5 illustrates the decadal trends for the different datasets over the 1979-1994 period. 
For the non-satellite datasets, the trends at 50 and 100 hPa are illustrated in panels (a) and (b), 
respectively; panel (c) illustrates the satellite-derived trends. The latitudes where the trends are 
statistically significant for the different datasets are listed in Table 5-4. [The Oort data, which 
has been used widely (e.g., Hansen et al., 1995; Santer et al., 1996), is not included in this plot 
owing to the fact that it spans a shorter period of time (1979-1989) than the other datasets.]. In 
the case of the MSU and Nash (SSU 15X) satellite data, the trend illustrated in panel (c) is in- 
dicative of a response function that spans a wide range in altitude (Figure 1); e.g., for MSU, 
about half of the signal originates from the upper troposphere at the low latitudes. Because of 
this, caution must be exercised in comparing the magnitudes of the non-satellite trends in pan- 
els (a) and (b) with those for the satellite in panel (c). This aspect could explain, in part, the 
lesser cooling obtained by the satellites relative to radiosondes in the tropical regions; however, 
this argument is contingent upon the trends in the tropical upper troposphere (not investigated 
in this report). The MSU indicates less cooling than Nash in the tropics. One reason for this 
could be that the Nash peak signal originates from a slightly higher altitude than MSU; again, 
though, the extent of the cooling/ warming trend in the upper troposphere needs to be consid- 
ered for a full explanation. The results are statistically insignificant in almost all of the datasets 
at the low latitudes. This could be in part due to the variable quality of the tropical data. It is 

conceivable that the radiosonde trends are significant over selected regions where the data are 

l-i tL 

reliable over long time periods, but that the significance aspect is destroyed whemreliable and 
unreliable data are combined to get a zonal mean. 

All datasets indicate a cooling of the entire northern hemisphere and the entire low and mid- 
latitude southern hemisphere at the 50 hPa level over this period. At the 100 hPa level, there is 
a cooling over most of the northern and southern latitudes. The midlatitude (30-60N) trends in 
the northern hemisphere exhibit a statistically significant (Table 5-4) cooling at both 50 and 100 
hPa levels, with the magnitude in this region being -0.5- IK/ decade. This feature is true for the 
satellite data as well. The similarity of the magnitude and significance in the mid-northern 
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hemisphere latitudes from the different datasets is particularly encouraging and suggests a ro- 
bust trend result for this time period. The trends in the southern hemisphere midlatitudes (-15- 
45S) range up to -0.5- IK/ decade, but are generally statistically insignificant over most of the 
area in almost all datasets, except Reanal. Note that the southern hemisphere radiosonde data 
has more uncertainties owing to fewer observing stations and data homogeneity problems (see 
5.2.5. 1). The non-satellite data indicates a warming at 50 hPa but a cooling at 100 hPa at the 
high southern latitudes, while the satellites indicate a cooling trend. Thus, as for the tropical 
trends, satellite-radiosonde intercomparisons in this region have to consider carefully the vari- 
ation of the trends with altitude (5.2.5. 3). The lack of statistically significant trends in the south- 
ern high latitudes need not imply that signifcant trends do not occur during particular seasons 
(e.g., Antarctic spring time). The high northern latitudes indicate a strong cooling (IK/ decade 
or more) in the 50 hPa, 100 hPa and satellite datasets. However, no trends are significant pole- 
ward of -70N owing to the large interannual variability there. There is a general consistency of 
the analyzed datasets (CPC, GSFC, Reanal) with trends derived directly from the instrumental 
data. Considering all datasets, the global lower stratospheric cooling trend over the 1979-1994 
period is estimated to be -0.6 K / decade. The 50-100 hPa cooling is consistent with earlier 
WMO results based on shorter records (1988 (e.g., Figure 6.17); 1990 (e.g.. Figure 2.4-5)) . 

Figure 6 shows the annual-mean trend over 1966-1994 at 50 hPa and comprises only the 
radiosonde record. Note that the Oort time series extends only through 1989. The cooling 
trends in the northern high latitudes, and in several other latitude belts, are less strong in the 
radiosonde datasets when the longer period is considered (note that the Berlin radiosonde time 
series for July, Figure 4, exhibits a similar feature). The cooling trend in the 30-60N belt is 
about 0.3K/ decade. The strong cooling trend in the Oort data in the high southern latitudes is 
consistent with Oort and Liu (1993), Parker et al. (1997) and the Angell data (D. Gaffen, per- 
sonal communication). Regions of statistically significant trends in the datasets are listed in Ta- 
ble 5-5. In the southern hemisphere, the two global radiosonde datasets indicate a significant 
cooling over broad belts in the low and midlatitudes, with the Oort data exhibiting this feature 
at even the higher latitudes. The Oort global-mean trend is -0.33 K/ decade over the 1966-1989 
period. In the northern hemisphere, again, the midlatitude regions stand out in terms of the sig- 
nificance of the estimated trends. Latitudes as low as 10-20 degrees exhibit significant trends 
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over the longer period considered. 


2.3.3 VERTICAL PROFILE 

Figure 7 shows the vertical and latitudinal structure of the zonal, annual-mean temperature 
trend as obtained from the SSU satellite measurements. The plot is constructed by considering 
5 km. thick levels from linear combinations of the weighting functions of the different SSU 
channels (Figure 1). Figure 7 shows panels with and without the inclusion of the volcanic pe- 
riods (i.e., the ‘no-volcano’ calculations omit two years of data following the El Chichon and 
Pinatubo eruptions; see Figure 2). The omission of the volcano-induced warming period (par- 
ticularly that due to Pinatubo near the end of the record) yields an enhanced cooling trend in 
the lower stratosphere. The vertical profile of the temperature trend in the middle and upper 
stratosphere between ~60N and -60S shows a strong cooling, particularly in the upper strato- 
sphere (up to 3K / decade). Cooling at these latitudes in major portions of the middle and upper 
stratosphere are seen to be statistically significant. 

We next focus on 45N latitude. At this latitude, lidar records from Haute Provence are avail- 
able, which afford a high vertical resolution above 30 km., relative to the other instrumental 
data available. Figure 8 displays the annual-mean trend profile for 1979-1998 updated from 
Keckhut et al (1995) and shows a cooling over the entire altitude range 35-70 km. of ~1-2K/ 
decade, but with statistical significance obtained only around 60 km. The vertical gradient in 
the profile of cooling between 40-50 km. differs somewhat from the 1979-1990 summer trend 
reported in WMO (1992, Figure 2-20). 

Figure 9 compares the vertical pattern of temperature trend at 45N obtained from different 
datasets for the 1979-1994 period. There is a broad agreement in the cooling at the lower strato- 
spheric altitudes, reiterating Figure 5. The vertical pattern of the trends from the various data 
are also in qualitative agreement, except for the lidar data (which, in any case, is not statistically 
significant over that height range). Note that the lidar trend for the 1979-1998 period (Figure 
8) exhibits better agreement with the satellite trend in Figure 9 (bottom panel), indicating a sen- 
sitivity of the decadal trend to the end year considered. Generally speaking, there is an approx- 
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imately uniform cooling of about 0.75 K/decade between ~80 and 5 hPa (-18-35 km.), 
followed by^ increasing cooling with height (e.g., -2.5 K/decade at 1 hPa {-50 km.}). The 
“analyses” datasets, examined here for p>10 hPa, are in approximate agreement with the in- 
strument-based data. 

2.3.4 ROCKET DATA AND TRENDS COMPARISONS 

Substantial portions of the three available rocket datasets (comprising US, Russian and Jap- 
anese rocketsondes; Table 5-2) have been either re-analyzed or updated since the review by 
Chanin (1993). Golitsyn et al. (1996) and Lysenko et al. (1997) updated the data from five dif- 
ferent locations over the period 1964-1990 or 1964-1995. Golitsyn et al. (1996) conclude a sta- 
tistically significant cooling from 25 to 75 km., except around 45 km. (Figure 10A). Lysenko 
et al. (1997) obtain a similar vertical profile of the trend for the individual rocketsonde sites 
(Figure 10B). They find a significant negative trend in the mesosphere particularly at the mid- 
latitude sites. 

Out of the 22 stations of the US network, 9 provide data for more than 20 years. Some series 
have noticeable gaps which prevents them from being used for trend determination. Indepen- 
dently two groups have revisited the data: Keckhut et al (1998) selected 6 low latitude sites (8S- 
28N), and Dunkerton et al (1998) selected 5 out of 6 of the low latitude sites plus Wallops Is- 
lands (37N). Both accounted for spurious jumps in the data by applying correction techniques. 

Keckhut et al. (1998) found a significant cooling between 1969 and 1993 of about 1-3 K / 
decade between 20 and 60 km. (Figure 1 1 A). A similar result is seen in the data available since 
1970 from the single and still operational Japanese rocket station (Figure 1 IB, updated from 
Komuro, 1989). Dunkerton et al. (1998), using data from 1962 to 1991, infer a downward trend 
of -1.7 K/decade for the altitude range 29-55 km.; they also obtain a solar-induced variation of 
-IK in amplitude. It should be noted that the amplitude of the lower-mesosphere cooling ob- 
served in the middle and high latitudes from the Russian rockets (between 3 and lOK/decade) 
is somewhat larger than from the US dataset. Nevertheless, all the four rocketsonde analyses 
shown in Figures 10 and 1 1, together with the lidar trend at 44N (Figure 8) are consistent in 
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yielding trends of -1-2K/ decade between about 30-50 km. 


Figure 12 compares the vertical profiles of trends over the 1979-1994 period from various 
datasets at ~30N, including the rocket stations at 28N (Cape Kennedy) and 34N (Point Mugu). 
As at 45N (Figure 9), almost all the datasets agree in the sign (though not in the precise mag- 
nitude) of temperature change below about 20 hPa (-27 km.). Above 10 hPa (-30 km.), both 
satellite and rocket trends yield increasing cooling with altitude, with a smaller value at the 28N 
rocket site. At 1 hPa, there is considerable divergence in the magnitudes of the two rocket 
trends. Because the rocket trends are denved from time series at individual locations, this may 
explain their greater uncertainty relative to the zonal-mean satellite trends. In a general sense, 
the vertical profile of the trend follows a pattern similar to that at 45N (Figure 9). 


2.4 Latitude-season trends 

The monthly, zonal-mean trends in the lower stratosphere are considered next. Figure 13 
illustrates the lower stratospheric temperature trends derived from MSU data over the period 
January 1979-May 1998. This represents an update from WMO (1995). Substantial negative 
trends are observed in the midlatitudes of both hemispheres during summer (- 0.5 K / decade), 
and in springtime polar regions (up to -3 to -4K/ decade). There is a broad domain of statistical 
significance in the cooling at the midlatitudes; in the NH, this occurs from -June through 
October while, in the SH, it occurs from about mid-October to April. Little or no cooling is 
observed in the tropics in Figure 13, although, as already stated, the MSU measurements in the 
tropics originate from a broad layer 50-150 hPa and thus implicitly include a trend signal from 
the upper troposphere (the weighting function maximum is near the tropical tropopause; Figure 
1). Thus, the tropical MSU result may not be indicative of a purely lower stratospheric trend 
(section 2.3.2; Figure 5). There is some symmetry evident in the observed trends in the two 
hemispheres, although most of the regions have values not significantly different from zero. 
Comparison with the MSU trends derived from data over 1979-1991 (Randel and Cobb, 1994; 
see Figure 1 1 in WMO, 1995) shows similar patterns in the SH, but substantial differences in 
the NH. The 1979-1991 NH trends peaked in winter midlatitudes, and the springtime polar 
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cooling was not statistically significant. Addition of data through mid- 1998 diminishes the 
northern midlatitude winter trends, and increases those over the pole in spring. Thus, owing to a 
high degree of interannual variability, especially during polar spring, the MSU trend estimates 
are dependent on the choice of the end year. 


The larger springtime polar trends are strongly influenced by the occurrence of several 
relatively cold years since about 1990, as shown in Figure 14 (top panel) for 80N using the 


NCEP re-analyzed data and MSU observations. Note that the NCEP reanalyses data have not /u*. ^ 
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been adjusted for the discontinuities due to switch-over to satellite temperature retrievals (1976' 1 ^ ^ 


for SH and 1979 for NH) which leads to a discontinuity in the analyzed temperatures for these 
years, associated with the model’s systematic error. Naujokat and Pawson (1996) have also 
noted the relatively cold 1994/1995 and 1995/1996 northern polar winters. The 1990s-averaged 
temperature in March is lower than that for the 1980s. Manney et al. (1996) suggest that an 
unusual circulation pattern (low planetary wave forcing and an intense polar vortex) was 
primarily responsible for this feature. In particular, March 1997 was very cold, about 18K below 
the 1980-1989 decadal average. Indeed, the monthly-mean 30 hPa temperature over the pole in 
March 1997 may have been one of the coldest years since 1966 (“Monthly Report on Climate 
System, 1997” by Japan Meteorological Agency). Newman et al. (1997) indicate that this 
coincided with the occurrence of the lowest total ozone amount in the northern hemisphere for 
this season. The extreme low temperatures in March 1997 appear to be associated with record 
low planetary wave activity (Coy et al., 1997). There are suggestions that the observed Arctic 
polar stratospheric conditions in recent years may be linked to changes in tropospheric 
circulation (s ee Cha pter 12). It is not certain whether these are secular trends or whether they are 
the consequence of a decadal-scale variability in the climate system. With regard to the 
frequency of major sudden stratospheric warmings, Labitzke and van Loon (1992) note that, 
before 1992, the largest lapse of time between two major warmings was about 4 years. In 
comparison, no major warming appears to have occurred between 1991 and 1997 i.e., over 
seven winters. This is broadly consistent with the sense of the temperature trend in the 1990s up 


to 1997. 


In the Antarctic, analysis of long-term records of radiosonde data continue to reveal sub- 
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stantial cooling trends in spring, further endorsing the early study of Newman and Randel 
(1988). This cooling trend is closely linked to springtime ozone depletion (Angell, 1986; Chu- 
bachi, 1986; Trenberth and Olsen, 1989; Jones and Shanklin, 1995; Butchart and Austin, 1996). 
This springtime cooling is an obviously large feature of the MSU temperature trends shown in 
Figure 13. Figure 14 (bottom) shows time series of 100 hPa temperatures in November over 
Halley Bay, derived from radiosonde data, together with time series from NCEP reanalysis in- 
terpolated to Halley Bay. Both these records reveal a decadal-scale change in temperature be- 
ginning in about the early 1980s, together with significant year-to-year variability (see also 
Chapter 4). The timing of the change in the Antarctic may be compared to that at 80N (March) 
around 1990, as illustrated in the top panel of Figure 14. It is notable that temperatures in NH 
springtime polar regions indicate a strong cooling in recent years and an enhancement above 
natural variability (also, compare Figure 13 with Figure 8-11 of WMO, 1995) that is reminis- 
cent of that observed in the Antarctic about a decade ago. 

The structure of the decadal-scale temperature change over Antarctica as a function of al- 
titude and month is shown in Figure 15, calculated as the difference in decade means (1987-96 
minus 1969-79), and averaged over seven radiosonde stations (Randel and Wu, 1998). These 
data show a significant cooling (of -6K) in the lower stratosphere in spring (October-Decem- 
ber). Significant cooling persists through austral summer (March), while no significant temper- 
ature changes are found during winter. These data also show a statistically significant warming 
trend (3K or more) at the uppermost data level (30 hPa; 24 km.) during spring. 

Kokin and Lysenko (1994) have analyzed the seasonal trends in the middle and upper 
stratosphere from their 5 rocketsonde station data for the 1972-1990 period. Generally, the 
cooling is evident almost throughout the year at all sites but there are exceptions. Figure 16 in- 
dicates a significant positive trend during springtime above Moldezhnaya (at -35 km.) and dur- 
ing winter above Volgograd and Balkash (at -40 km.). The warming feature for Moldezhnaya 
is consistent with the springtime warming inferred for slightly lower altitudes in the Antarctic 
by Randel (1988), and with that shown in Figure 15. The warming is located above the domain 
of cooling caused by the lower stratospheric ozone depletion. It is interesting that such an effect 
also takes place during the northern winter over Volgograd and Balkash locations. 
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2.5 Uncertainties in trends estimated from observations 


Determining stratospheric temperature trends from long-term observations is complicated by 
the presence of additional, non-trend variability in the data. Two types of phenomena contribute to 
the uncertainty in trend estimates. The first is true atmospheric variability that is not trend-like in 
nature. Major sources of such variability include: the (quasi-) periodic signals associated with the 
annual cycle, the quasi-biennial oscillation, the solar cycle, and the El Nino-Southern Oscillation 
(ENSO). In addition, stratospheric temperatures vary in response to episodic injections of volcanic 
aerosols. To first approximation, these atmospheric phenomena have negliglible effects on the long- 
term temperature trend because they are periodic or of relatively short duration. Nevertheless, 
because current data records are only a few decades long, at most, these phenomena may appear to 
enhance or reduce an underlying trend. At a minimum, the additional temperature variability 
associated with these signals reduces the statistical confidence with which long-term trends can be 
identified. While periodic signals can be removed, the effects of sporadic events are more difficult 
to model and remove. Furthermore, there may be long-term trends in these cycles and forcings that 
confound the analysis. A second source of uncertainty is due to spurious signals in the time series 
that are due to changes in methods of observation rather than to changes in the atmosphere. The 
problem of detecting temperature trends in the presence of changes in the bias characteristics of v 
the observations is receiving increased attention (Christy, 1995; Santer et al., 1998). It seems likely 
that over the next few years better methods will be employed to quantify and reduce the uncertainty 
in stratospheric temperature trend estimates attributable to these spurious signals. 

2.5.1 UNCERTAINTIES ASSOCIATED WITH RADIOSONDE DATA 


Although most radiosonde analyses show cooling of the lower stratosphere in recent decades, 
it is important to recognize that they all rely on subsets of the same basic dataset, the global 
observing system upper-air network. This network is fundamentally a meteorological one, not a 
climate monitoring network, not a reference network for satellite observations, and not a network 
for detection of stratospheric change. Whatever difficulties plague the radiosonde network when it 
is construed for temperature trends analysis will affect all analyses that use those data. 
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Karl et al. (1995) and Christy ( 1 995) have reviewed some of the problems associated with using 
radiosonde data for the detection of atmospheric temperature trends. These fall into two categories: 
the uneven spatial distribution of the observations, and temporal discontinuities in station records. 

The radiosonde network is predominantly a northern hemisphere, midlatitude, land network. 
About half the stations are in the 30-60N latitude band, and less than 20% are in the southern 
hemisphere (Oort and Liu, 1993). Moreover, the uneven distribution of stations is worse for 
stratospheric data than for the lower troposphere, because low latitude and southern hemisphere 
soundings have a higher probability of taking only one observation daily (other stations make two, 
and many formerly made four), and because the soundings more often terminate at lower altitudes 
(Oort and Liu, 1 993). Estimates of layer-mean trends, and comparisons of trends at different levels, 
are less meaningful when data at the top of the layer are fewer than at the bottom. 

For trend detection, the temporal homogeneity of the data is the key. As discussed by Parker 
(1985), Gaffen (1994), Finger et al. (1993), and Parker and Cox (1995), numerous changes in 
operational methods have led to discontinuities in the bias characteristics of upper-air temperature 
observations, which are particularly severe in the lower stratosphere (Gaffen, 1994). Because 
radiosondes are essentially expendable probes (although some are recovered and reconditioned), 
changing methods is much easier than, for example, in surface observations at fixed locations. 
Effects on data have been shown for changes in instrument manufacturer and replacement of old 
models with newer ones from the same manufacturer (Parker, 1985; Gaffen, 1994), changes in time 
of observation (Elliott and Gaffen, 1991; Zhai and Eskridge, 1996), changes in the lag 
characteristics of the temperature sensors (Parker, 1985; Huovila and Tuominen, 1990), and even 
changes in the length of the suspension cord connecting the radiosonde balloon with the instrument 
package (Suzuki and Asahi, 1978; Gaffen, 1994) and balloon type (Parker and Cox, 1995). 

Daytime stratosphenc temperature data in the early years of radiosonde operations were 
particularly affected by errors due to solar radiation, and substantial changes in both data correction 
methods (e.g., Sera se, 1956; Teweles and Finger, 1960) and instrument design have been made to 
address the problem. In general, the result has been a reduction in a high bias over time, leading to 
an artificial “cooling” in the data (Gaffen, 1994). The example in Figure 17, showing monthly 
mean 100 hPa temperatures at Tahiti, shows that a 1976 change from one model of the French 
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Mesural radiosonde to another, each with different temperature sensors, was associated with an 
artificial temperature drop of several degrees. Given estimated temperature trends on the order of 
tenths of a degree K per decade, such inhomogeneities introduce substantial uncertainty regarding 
the magnitude of the trends in the lower stratosphere. Furthermore, Luers (1990) has demonstrated 
that daytime radiosonde temperature errors can exceed IK at altitudes above 20 km (50 hPa) and 
that the magnitude of the error is a strong function of the temperature and radiative environment, 
which suggests that, as the atmosphere changes, so will the nature of the measurement errors. 

A few investigators have attempted to adjust radiosonde temperature time series to account for 
“change-points,” or level shifts like the one illustrated in Figure 17. Miller et al. (1992), using data 
from 1970-1986, made adjustments to lower stratospheric data at four of 62 stations in Angell’s 
(1988) network on the basis of a statistical regression model that includes a level shift term. 

However, Gaffen ( 1 994) concluded that, over a longer period, many more of the Angell stations 
showed data inhomogeneities. Parker et al. (1997) have made adjustments to temperature data from 
Australia and New Zealand for the period 1979-1995 by using MSU data as a reference time series, 
and station histories (Gaffen, 1 994) to identify potential change-points. The adjustments (of earlier 
data relative to 1995 data) ranged from 0 to -3.3 K and reduced the estimated zonal-mean 
temperature change between the periods 1987-1996 and 1965-1974 at about 30S (30 hPa) from -2.5 
K/ decade to about -1.25 YJ decade (Parker et al., 1997). 

2.5.2 UNCERTAINTIES ASSOCIATED WITH ANALYSES OF SSU SATELLITE 

DATA A/ — i/i j 

The stratospheric temperature analyses from CPC (NCEP) are an operational product, 
derived using a combination of satellite and radiosonde temperature measurements. Radiosonde 
data contribute to the analyses in the NH over the 70-10 hPa levels; satellite data alone are used 
in the tropics and SH, and over the entire globe above 10 hPa (Gelman et al., 1986; section 2.1). 

The satellite temperature retrievals are from the TIROS Operational Vertical Sounder (TOVS), 
which have been operational since late 1978. A series of TOVS instruments have been put into 
orbit aboard a succession of operational satellites; these instruments do not yield identical 
radiance measurements for a variety of reasons, and derived temperatures may change 
substantially when a new instrument is introduced (Nash and Forrester, 1986). Finger et al. 

(1993) have compared the operationally derived temperatures with co-located rocketsonde and 
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lidar observations, and find systematic biases of order +/- 3 to 6 K in the upper stratosphere (due 
primarily to the low vertical resolution of TOVS). These biases furthermore change with the 
introduction of new operational satellites, and Finger et al (1993) provide a set of recommended 
corrections to the temperature data (dependent on the particular satellite instrument for each 
time period) which have been used by CPC. 

In spite of the application of the adjustments recommended by Finger et al. (1993), time 

series of temperature anomalies from the CPC analyses still exhibit significant discontinuity 

near the times of satellite transitions. Figure 18 shows deseasonalized CPC temperature 

anomalies over the equator at 1 hPa, together with a time series of equatorial anomalies from the 

overlap-adjusted (i.e., accounting for the discontinuity in satellites) SSU Channel 27 data 

(whose weighting function peaks near 1.5 hPa; Figure 1). These time series show very different 

characteristics, with apparent discontinuity or 'jumps' in the CPC data which are coincident with 

satellite changes. Two particularly large changes are seen in Figure 18 in mid- 1984 and late 

1988; the specific nature of these discontinuities depends on latitude and altitude. Overall, the 

presence of such discontinuities in the CPC data limits their reliability for robust estimates of 

upper stratospheric trends (as in, say. Figure 9) and is even problematic for the determination of 

interannual variability .^The overlap-adjusted SSU data (see F igur e 2) used in Figur es 7, 9 and 12 

do not exhibit such obvious problems. - wc vAv "lA /Vt. 1 1 * ( $ 
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2.5.3 UNCERTAINTIES IN SATELLITE-RADIOSONDE TREND 
INTERCOMPARISONS 

A distinct advantage of the satellite instruments over in situ ones is their globally extensive cov- 
erage. However, as already mentioned, this is tempered by the fact that the signals that they receive 
originate from a broad range of altitudes (Figure 1). This is in contrast to the specific altitudes of 
measurements in the case of the ground-based instruments located at specific sites. This feature of 
the satellite trends complicates the interpretation for any particular vertical region of the atmo- 
sphere and, more particularly hampers a rigorous comparison with, say, the radiosonde trends (see 
section 2.3.2). As an example, consider the problem of the lower stratospheric temperature trends. 
The MSU’s Channel 4 ‘senses’ the entire extent of the lower part of the stratosphere, and even the 
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upper troposphere at low latitudes. This poses problems in the precise intercomparison of presently 
available satellite-based trends with those from ground-based instruments. In the tropics, approxi- 
mately half of the signal originates from the upper troposphere, leading to a potential misinterpre- 
tation of the actual lower stratospheric temperature trend based on MSU alone. This problem can 
become acute particularly if the tropical upper troposphere and lower stratosphere have tempera- 
ture trends of opposite signs. A similar comment also applies to the interpretation of the strato- 
sphere trends from SSU measurements (Figure 1). Further, because of the areal coverage of the low 
latitudes, the global-means from satellite data and those from the m-situ instruments may be com- 
parable only after appropriate adjustments are made for the differential sampling by the two kinds 
of instruments. Therefore, caution must be exercised in the interpretation of satellite-based trends 
vis-a-vis radiosonde and other ground-based instruments. Besides, satellite data interpretations 
also have to cope with problems involving temporal discontinuity, instrument calibration and orbit 
drift. 


2.5.4 UNCERTAINTIES ASSOCIATED WITH ROCKET DATA 

The rocket data are very useful as they were the only observations of the 30-80 km. region 
before the lidars started operating. However, determining quantitative trends from rocket data is 
complicated by both physical and measurement issues. A first difficulty with the rocket data is 
that there have been instrumental changes and the measurements come from different types of 
sensors (Arcasonde, datasonde, falling spheres). However, Dunkerton et al (1998) have found 
that these changes were a less important source of error than previously suggested. The major 
source of error, and the origin of the observed spurious jumps, seems to be due to the change of 
corrections of the data to take into account aerodynamic heating. Most of the earlier analyses did 
not take full account of the changes and of the spurious jumps in the data that ensued from the 
above-mentioned difficulties. These points have been considered in-depth by Keckhut et al. 
(1998) and Dunkerton et al. (1998), which resulted in a very limited number of US stations that 
could be used for determining trends. Yet another source of uncertainty is due to the different 
time of measurements, as the amplitude of tidal influence may not be negligible at these 
altitudes (+/- 2K around 40-45 km. according to Gille et al. (1991) and Keckhut et al. (1996)). 
This may explain the small error bar in the Ryori Japanese measurements (Figure 1 IB), which 
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are always conducted at the same local time. The factor related to the local measurement time 
was accounted for in the analysis of the US rocket data by Keckhut et al. (1998). Since rocket 
data are available only from a few locations around the globe, there is a difficulty in ascertaining 
consistency of the trend and its significance when compared against zonal-mean satellite data. 

2.5.5 UNCERTAINTIES ASSOCIATED WITH LIDAR RECORD 

Temperature obtained from Rayleigh lidar are affected by the presence of aerosols. Thus, a 
measure of the atmospheric scattering ratio is required after major volcanic eruptions. 

Otherwise, the temperatures are given in absolute value as a function of altitude from 30 km. 
upwards, without any need of external calibration. An accuracy of 1% is easily attained, with a 
principal limit for ascertaining the significance of a trend being the length of the available data 
set. Using the actual measurements at the Haute Provence site, it was found that the 
establishment of a significance in the trend at the 95% level in the upper stratosphere required 
20.5 years of data for summer and 35 years for winter trends. More years are required for the 
wintertime owing to the increased variability present in that season (Keckhut et al. 1995). Of 
course, the length of a period needed to establish statistical significance also depends on the 
amplitude of the signal. Compared to rocket data which are made at a specific time (even though 
it may not be the same local time for different sites), the lidar data can be made at any time 
during the night. This could constitute a potential -source of uncertainty owing to tidal effects, 
and may explain, in part, the large error bars observed in Figure 8. This could be improved by 
selecting data corresponding to the same local time. 


2.6 Issues concerning variability 

2.6.1 VOLCANIC AEROSOL INFLUENCES 

There are a number of factors that influence temperature trends in the stratosphere. Among 
the most significant as a short-term phenomenon is the large build-up of stratospheric aerosols 
following volcanic eruptions (see Chapter 3). The SAGE extinction and other satellite data 


25 



(e.g., McCormick et al., 1995) reveal a notable increase in stratospheric aerosol concentrations 
for 1-2 years following volcanic eruptions. Over the past two decades, such transient enhance- 
ments have come about due to volcanic eruptions of different intensities. Data from earlier 
ground-based observations reveal other episodes of volcanic loading of the stratosphere dating 
back to the 1960s. It is well understood that aerosols injected into the lower stratosphere by ma- 
jor volcanic eruptions result in a warming of this region of the atmosphere owing to enhanced 
absorption of solar radiation and the upwelling terrestrial infrared radiation (Pollack and Ack- 
erman, 1983; WMO, 1988; Labitzke and McCormick, 1992). 

Three volcanos stand out in particular since the 1960s, the period when widespread and 
routine radiosonde observations began. The three volcanos of particular importance for climate 
variations have been - Agung (Bali, Indonesia, 1963), El Chichon (Mexico, 1982) and Pinatubo 
(Philippines, 1991). In each one of these three cases, well-documented instrumental records in- 
dicate that the temperatures increased and stayed elevated till the aerosol concentrations were 
depleted completely (WMO: 1988, 1995). All three eruptions produced somewhat similar 
warming characteristics (Angell, 1993) viz., a warming of the lower stratosphere (-15-25 km.) 
centered in the tropics (-30N-30S), with a magnitude of 1-2K (up to 3K for Pinatubo). The 
warming diminishes in time, but anomalies are observed for approximately 2 years following 
each eruption (Figures 2 (SSU 1 5X) and 3). Because the volcanic effects are episodic and clear- 
ly identifiable, their direct effect on calculation of trends is minimal (the volcanic time periods 
simply need to be omitted prior to trend calculation), except as they possibly overlap with other 
potentially causal phenomena such as the solar cycle. 

The magnitude and evolution of the transient warming agree reasonably well between the ra- 
diosonde and satellite datasets (Labitzke and McCormick, 1992; Christy and Drouilhet, 1994; Ran- 
del and Cobb, 1994; Figure 3 (top) and (bottom) panels). Relative to the analyses in WMO (1995), 
the temperature of the global-mean lower stratosphere has become progressively colder in both the 
radiosonde and satellite time series (Figures 2 (bottom panel) and 3). 

In the upper stratosphere and mesospheric regions, an indirect effect could be expected from 
the change of circulation and J or the upwelling flux arising from the lower stratospheric heating. 
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Such an indirect effect appears to have been observed after the Pinatubo eruption in the OHP lidar 
data. A cooling of 1 .5K in the upper stratosphere and a warming of 5K at 60-80 km. were detected 
(Keckhut et al., 1995). Thus, a trend analysis for these regions also requires the omission or cor- 
rection of post- volcanic eruption data (e.g., such a correction was applied to the data plotted in Fig- 
ure 8 by adding a term proportional to the aerosol optical thickness in the analysis technique). 

2.6.2 SOLAR CYCLE 

The 1 1-year modulation of the UV solar flux which is now well documented is expected to, 
through photochemistry, influence stratospheric ozone and therefore stratospheric temperature. 

A number of investigators have attempted to identify a signature of the 1 1-year solar cycle in the 
temperature dataset. The solar proxy used in most studies is the 10.7 cm. radio flux which spans 
the longest period, even though more realistic proxies would be the He I line, Mg II line or UV 
irradiance. However, several analyses using these different proxies justify the choice of the 10.7 
cm. flux (Donnelly et al., 1986; Keckhut et al., 1995). 

Satellite data provide a global view of the signature due to solar variations, but the time 
series is relatively short. The overlap-adjusted SSU and MSU data sets (spanning 1979-1995) 
exhibit a coherent temperature variation approximately in phase with the solar cycle. Figure 19 
shows the vertical-latitudinal structure of the solar signal (derived via regression onto the 10.7 
cm. solar flux time series) using the overlap-adjusted MSU and SSU data. Even though the 
dataset is limited to 17 years, this shows a statistically significant solar component of order 0.5 - 
1 .0 K throughout most of the low latitude (30N-30S) stratosphere, with a maximum near 40 km. 
The spatial patterns show maxima m the tropics, with an approximate symmetry about the 
equator. A small solar response is observed in these data at the high latitudes. The 0.5- IK solar 
signal seen in the Nash data is in reasonable agreement with results from the longer records of 
radiosondes/ rocketsondes discussed below. At northern midlatitudes, the satellite-derived 
signature, which is not statistically significant, goes from slightly positive at -17 km. to slightly 
negative at -25 km. with a null value at 45-50 km. This is similar to the OHP lidar record 
(Keckhut et al., 1995), and to the Volgograd rocket record (Kokin et al., 1990). At high 
latitudes, although the satellite data uncertainty is large, there is a hint of a large positive 
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response in the mesosphere, as is observed in the Heiss Island rocket data (Kokin et al, 1990). 


The longest time series is provided by radiosondes and rockets; even then, they barely cover 
at best 4 solar cycles. There is, thus, in general, some difficulty in establishing a firm 
relationship with the 1 1-year solar cycle. Labitzke and van Loon (1989) and van Loon and 
Labitzke (1990) use the Berlin meteorological analyses, beginning in 1964 and spanning the NH 
lower stratosphere, to isolate coherent temperature cycles in the subtropics in phase with the 
solar cycle. Labitzke and van Loon (1997) update their previous analysis and find a correlation 
of 0.7 with the 30 hPa geopotential height (which can be taken to imply a similar correlation 
with layer-mean temperature below this level). The maximum correlation occurs over W. 

Pacific near China. Using the NCEP reanalyzed global dataset, van Loon and Labitzke (1998) 
show a similar solar signal manifest in the Southern Hemisphere, with higher correlations in the 
tropical regions. 

Angell (1991b) has used radiosonde and rocketsonde data to deduce tropical and NH 
midlatitude solar cycle variations of approximately 0.4 and 0.8 K in the lower and upper 
stratosphere, respectively (changes per 100 units of 10.7 cm. solar flux, or approximately solar 
maximum minus solar minimum values). Dunkerton and Baldwin (1992) isolate a weak solar 
cycle in CPC temperatures in the NH winter lower stratosphere, using data from 1964-1991 
(these analyses are based primarily on radiosonde data). Isolation of a solar cycle signal in CPC 
upper stratospheric temperature data beginning in 1979 is somewhat problematic, due to the 
spurious discontinuities introduced by satellite changes (see Fig. 18; section 2.5.2). This 
suggests that the solar cycle variations derived from these data (e.g. Kodera and Yamazaki, 

1990; Hood et al., 1993) should be treated with caution. 

In their rocket data analysis, Dunkerton et al. (1998) found a 1.1K response to the solar cycle 
for the integrated altitude range of 29-55 km. (Figure 20). Kokin et al. (1990), Angell (1991b) 
and Mohanakumar (1995), using rocket data, and Keckhut et al. (1995), using lidar data, infer a 
clear solar signature in the mesosphere of +4 to 10K per 100 units of 10.7 cm. flux. On the other 
hand, the results obtained in the upper stratosphere and around the stratopause are different in 
both amplitude and sign for the different sites, and are also variable with season. 
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From the records, it thus appears that the solar cycle signature in stratospheric temperatures 
need not be uniform and identical all over the globe and at all altitudes. This has been manifest 
in Labitzke and van Loon (1997) on the horizontal scale, as well as on the vertical scale in 
Chanin and Keckhut (1991 ), and may be attributable to the role of planetary waves. An 
additional point to note is that the solar-induced temperature changes need not occur at the same 
latitudes as any changes in ozone. 

In the lower stratosphere, the effect of solar variations is found to have a relatively small 
effect on trend calculations for time series longer than 15 years as those analyzed here. 
Regression estimates of trends neglecting a solar cycle term change by only -10% (this 
sensitivity was determined by testing several of the time series analyzed here). Elsewhere in the 
stratosphere, the amplitude of the solar cycle signature has the potential to introduce a bias in 
trend estimates (especially if the number of cycles involved in the time series is small). The 
determination of the amplitude and structure of the solar cycle effect may be disturbed in the 
two recent cycles by the presence of two major volcanic eruptions spaced about nine years apart, 
which is not different enough from 1 1 years, and so has the potential to be confused with the 
solar cycle. 

2.6.3 QBO TEMPERATURE VARIATIONS 

Interannual variability of zonal winds and temperature in the tropical stratosphere is 
dominated by an approximate two year periodicity - the quasi-biennial oscillation (QBO). The 
QBO is most often characterized by downward propagating zonal wind reversals in the tropics 
(e.g. Naujokat, 1986; Reid, 1994), but the QBO is also prevalent in tropical temperatures (which 
are in thermal wind balance with the zonal winds). Aspects of the QBO in global stratospheric 
temperatures have been derived from radiosonde observations (Dunkerton and Delisi, 1985; 
Labitzke and van Loon, 1988; Angell, 1997) and satellite-derived data sets (Lait et al., 1989; 
Randel and Cobb, 1994; Christy and Drouilhet, 1994). Randel et al (1998) have used UKMO 
assimilated wind and temperature fields (Swinbank and O’Neill, 1994), together with overlap- 
adjusted SSU data, to isolate global patterns of the QBO. The results of these studies produce 


29 



the following picture of the QBO in temperature: 


1) Tropical QBO temperature anomalies maximize in the lower stratosphere (-20-27 km.) over 
10N-10S, with an amplitude of +/- 2 to 4 K. There is a secondary maximum in the tropical upper 
stratosphere (-35-45 km.) of similar magnitude, approximately out of phase with the lower 
stratosphere maxima. The structure in temperature is qualitatively similar to the 2-cell structure 
observed in ozone (see also Chapter 4 of WMO, 1998). 

2) There is a substantial midlatitude component of the QBO in temperature, maximizing over ~ 
20-50 degrees latitude, and approximately out of phase with the tropical anomalies. The 
midlatitude patterns also exhibit maxima in the lower and upper stratosphere, which are 
vertically out of phase. An intriguing aspect of the midlatitude QBO is that it is seasonally 
synchronized, occurring only during winter-spring of the respective hemispheres; this produces 
a highly asymmetric temperature pattern at any given time. 

3) There is also evidence for a coupling of the QBO with the lower stratosphere polar vortex in 
late winter-spring of each hemisphere (Holton and Tan, 1980, 1982; Lait et al., 1989). The polar 
QBO temperatures are out of phase with those in the tropics. Observational evidence for the 
polar QBO is less statistically significant than that in the tropics or midlatitudes, at least partly 
because of the high level of 'natural' interannual variability in polar vortex structure. Evidence of 
a quasi-biennial modulation of the southern hemisphere stratospheric polar vortex has been 
presented by Baldwin and Dunkerton (1998b). An interesting question arising from recent 
studies (Salby et al., 1997; Baldwin and Dunkerton, 1998a) is whether the interactions between 
the QBO and a so-called biennial oscillation, could be a factor in the observed signature of the 
solar cycle. 

Because the available satellite and radiosonde data records span many observed QBO cycles, 
it is possible to empirically model and isolate QBO variations in temperature with a good degree 
of confidence. The QBO has little effect on calculated trends, although trend significance levels 
are increased by accurate characterization of QBO variability. 


30 



2.6.4 PLANETARY WAVE EFFECTS 


The trend, solar cycle, volcanic and QBO variations account for the majority of observed 
interannual temperature variability in the tropics and low latitudes. Over middle and high 
latitudes there is a substantial amount of residual interannual variance, particularly during winter 
(e.g. Dunkerton and Baldwin, 1992). This additional variance is due to 'natural' meteorological 
fluctuations, in particular the presence or absence of stratospheric warming (or planetary wave) 
events during a particular month (a large wave event corresponds to a warm month, and vice 
versa). This relatively high level of natural variability of the polar vortex means that somewhat 
larger natural or anthropogenic signals (e.g. trends due to changes in radiative species, solar or 
QBO variations) are required for statistical significance. 

2.6.5 ENSO 

El Nino events are known to have a major impact on lower stratospheric temperature anom- 
alies in the tropics, with a warming of the sea surface temperature (SST) and the increase in 
moist convection being accompanied by a cooling of the lower stratosphere (Pan and Oort, 

1983; Reid, 1994) This inverse relationship between tropical surface and lower stratosphere 
temperatures is also manifest in the long-term record (Sun and Oort, 1995). Since tropical SSTs 
are known to have been increasing in recent decades (e.g. Graham, 1995), this suggests a po- 
tential contribution towards a tropical lower stratospheric cooling trend. It remains to be deter- 
mined how the secular trends in the upper troposphere and lower stratosphere are quantitatively 
related to SST variations and trends, both in the tropics and extratropics. 


2. 7 Changes in TVopopause Height 

Temperature changes in the upper troposphere and lower stratosphere may induce a change 
in the tropopause height, which may complicate the determination of temperature trends in the 
tropopause region. This has been first suggested by Fortuin and Kelder (1996) in the context of 
radiative cooling caused by the trends in ozone and other greenhouse gases. To date, most mod- 
els would be unable to represent such change as they are much smaller than the typical vertical 
spacing in models. Therefore, inferences have to rely entirely on the review of the available ob- 
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servations. Several results obtained at northern midlatitudes provide some indication of such 
changes, even though there is a lack of agreement on the amplitude of the effect. 

An increase of the tropopause height by 150 +/- 70 m (2-sigma) has been found at Hohen- 
peissenberg over the last 30 years (Steinbrecht et al., 1998). This increase is correlated with the 
decrease in the ozone mixing ratio in the lower stratosphere and with a tropospheric warming 
at 5 km. of 0.7 +/- 0.3 K (2-sigma) since 1967 (Figure 21). Bojkov and Fioletov (1997) indicate 
that the tropopause over Central Europe has increased by about 300 m over the last twenty 
years. For the tropics, Reid (personal communication) has carried out a preliminary analysis of 
data from Truk (7.5N, 151.8E) for the period 1965-1994, and has found no detectable change 
in the temperature but a significant change in height, consistent with a warming of the tropo- 
sphere (7.7 m/yr or 230 m in 30 years). These results have yet to be extended to include other 
stations and have not yet accounted for ENSO and volcanic influences. On a global scale, Hoin- 
ka ( 1 998a) has made a study of the tropopause surfaces based on European Centre for Medium- 
range Weather Forecasting reanalyses data, but for a shorter time period. His results indicate 
that the sign of the trend could vary depending on the latitude and longitude. A more recent 
work (Hoinka, 1 998b) analyzes the trends for tropopause height, potential temperature and wa- 
ter vapor mixing ratio. Clearly, this subject warrants more scrutiny on all spatial scales in order 
to quantify more rigorously the increase m the height of the tropopause, particularly over the 
past decade. 


3. MODEL SIMULATIONS 


3.1 Background 
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In this section, we discuss results from model investigations that have analyzed the effects 
due to various natural and anthropogenic factors, notably changes in trace species, upon strato- 
spheric temperature trends and variations. We focus m particular on the changes in the lower 
stratospheric region, and also discuss the vertical profile of the modeled trends from the lower- 
to-upper stratosphere. Where possible and relevant, available model results are compared with 
observations. 

Numerical models based on fundamental understanding of radiative, dynamical and chem- 
ical processes constitute essential tools for understanding the effects of specific mechanisms on 
temperature trends and variability in the stratosphere, and for interpreting observed tempera- 
ture changes in terms of specific mechanisms. The numerical models used thus far have at- 
tempted to include, to varying degrees, the relevant components of the climate system which 
could influence stratospheric temperatures. This includes interactions among radiative, chem- 
ical and dynamical processes The models also attempt to capture the important links between 
the stratosphere, troposphere and mesosphere. 

It is well-recognized that the global, annual-mean thermal profile in the stratosphere repre- 
sents a balance between solar radiative heating and longwave radiative cooling, involving 
mainly ozone, carbon dioxide, water vapor, methane, nitrous oxide, halocarbons and aerosols 
(see Goody and Yung, Chapter 9 and references therein). In the context of the general global 
stratosphere, dynamical effects also become a factor in determining the thermal profile. Since 
the late 1950s, with increasing knowledge of trace species’ concentration changes and their op- 
tical properties, numerical models have played a significant role in highlighting the potential 
roles of various constituents and the different mechanisms operating in the stratosphere. For ex- 
ample, WMO (1986; 1988) concluded that changes in the concentrations of trace gases and 
aerosols could substantially perturb the radiative balance of the contemporary stratosphere and 
thereby affect its thermal state. 

Early numerical models were developed as one-dimensional ones on the basis that the glo- 
bal, annual-mean stratosphere is in radiative equilibrium. Together with the assumption of a ra- 
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diative-convective equilibrium in the troposphere, this led to the so-called one-dimensional 
radiative-convective models (ID RCMs) which have been widely employed to study effects 
due to trace gas perturbations (WMO, 1986). The radiative-convective models represent the at- 
mosphere as a single global-mean vertical column (Manabe and Wetherald, 1967; Ramanathan, 
1981). The temperatures at the surface and in the troposphere are determined by both radiation 
and by some representation of the processes that determine the vertical advection of heat by the 
fluid motions (normally via a convection scheme). In the stratosphere, the thermal state deter- 
mined is a balance between absorbed solar radiation, and absorbed and emitted thermal infra- 
red radiation. Radiative-convective models can also include chemical processes to represent 
feedbacks between temperature and chemical constituents (e.g. Bruhl and Crutzen, 1988). 

A variation of the RCMs is the so-called Fixed Dynamical Heating model or FDH (Fels and 
Kaplan, 1975; Ramanathan and Dickinson, 1979; Fels et al., 1980). The contemporary FDH 
models (e.g., WMO, 1992) hold the tropospheric temperature, humidity and cloud fields fixed 
and allow for changes in the stratospheric temperature in response to changes in radiatively 
active species. It is assumed that, in the unperturbed state, the radiative heating is exactly 
balanced by the dynamical heating at each height. If the concentration of a radiatively-active 
constituent is altered then the radiative heating field is altered. It is assumed that the dynamical 
heating remains unchanged, and that the temperature field adjusts in response to the 
perturbation. In turn, this alters the radiative heating field such that it again exactly balances the 
dynamical field. The timescale for adjustment to a specific perturbation typically varies on the 
order of a few months in the lower stratosphere to a few days in the upper stratosphere. 

The application of the RCM and FDH model concepts for understanding stratospheric 
temperature changes has evolved with time (see WMO: 1990, 1992, 1995). Both types of 
models have been extensively used for gaining perspectives into the thermal effects due to the 
observed and projected changes in radiatively-active trace gases. These simple models, though, 
have important limitations. In particular, the FDH models can only predict temperature changes 
at any location due to constituent changes occurring there; i.e. the response is entirely localized 
within that particular stratospheric column. As will be seen later, changes in the circulation as a 
consequence of a radiative perturbation at any latitude can influence temperature changes at 
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other latitudes. Hence, it should not be expected that FDH or FDH-like models can realistically 
simulate features such as the latitudinal and monthly variation of trends. Nonetheless, these 
simple models remain useful in revealing reasonable, first-order solutions of the problem. 


There has been a steady progression from the simple RCMs and FDH models to the three- 
dimensional general circulation model (GCM; see WMO, 1986 for an early discussion of 
models used for studying the stratosphere) which seeks to represent the radiative-dynamical 
(and even chemical in some instances) interactions in their entirety. Such models have 
representations of radiative processes that may be less complete and accurate than in the ID 
models, but provide the best means of mimicking the real atmosphere. It may be noted that one 
useful aspect of the FDH simulations is that their comparison vis-a-vis a GCM simulation 
enables a quantitative distinction to be made between the radiative and radiative-dynamical 
responses of the stratosphere to a specific perturbation. 


3.2 Well-mixed greenhouse gases 

While increases in well-mixed greenhouse gases (carbon dioxide, methane, nitrous oxide, 
halocarbons) warm the surface, their effects on the lower stratospheric temperature varies, 
primarily because of the location of the absorption bands of each of the species and the radiative 
characteristics of the atmosphere in those bands. The essential radiative processes are that, first, 
an increase in greenhouse gas concentration enhances the thermal infrared emissivity of a layer 
in the stratosphere; hence, if the radiation absorbed by this layer remains fixed and other factors 
remain the same, then, to achieve equilibrium, the same amount of energy has to be emitted at a 
lower temperature and the layer cools. Second, an increase in emissivity due to increased 
concentration also implies an increase in the thermal infrared absorptivity; a larger amount of 
the radiation emitted by the troposphere will be absorbed in the stratosphere, leading to a 
warming tendency. Note that the upwelling/ downwelling radiation reaching the lower 
stratosphere is itself dependent on the infrared opacity and the thermal state of the troposphere/ 
middle-and-upper stratosphere. The net result is a balance involving these processes 
(Ramanathan et al., 1985; Ramaswamy and Bowen, 1994; Pinnock et al., 1995). 
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For carbon dioxide, the main 15 pm band is saturated over quite short distances. Hence the 
upwelling radiation reaching the lower stratosphere originates from the cold upper troposphere. 
When the C0 2 concentration is increased, the increase in absorbed radiation is quite small and 
the effect of the increased emission dominates, leading to a cooling at all heights in the 
stratosphere. For gases such as the CFCs, their absorption bands are generally in the 8-13 pm 
‘atmospheric window’ (WMO, 1986) - hence much of the upwelling radiation is originating 
from the warm lower troposphere. In this case the increase in the absorbed energy is generally 
greater than the increase in emitted energy and a warming of the lower stratosphere results 
(although there are exceptions to this case if the halocarbon absorption is in a region of already 
large atmospheric absorption - see Pinnock et al., 1995). Methane and nitrous oxide are mid-way 
between these two regimes so that an increase in concentration has little effect on the lower 
stratospheric temperature as the increased absorption and increased emission almost balance. In 
the upper stratosphere, increases in all well-mixed gases leads to a cooling as the increased 
emission effect becomes greater than that due to the increased absorption. 

One significant consequence of this differing behavior concerns the common use of 
‘equivalent C0 2 ’ in climate simulations. Equivalent C0 2 is the amount of C0 2 used in a model 
calculation that results in the same radiative forcing of the surface-troposphere system as a 
mixture of greenhouse gases (see e.g., IPCC, 1995). Although the equivalence concept works 
reasonably well for tropospheric climate change (IPCC, 1995), it does not work well for 
stratospheric temperature changes because increases in different greenhouse gases cause 
different stratospheric temperature responses, thus requiring the consideration of each gas 
explicitly. In fact, some trace gases warm, notcool like C0 2 , the lower stratospheric region near 
the tropopause (WMO, 1986). Even if all gases were to cause an effect of the same sign in the 
stratosphere, an equivalence with respect to C0 2 that holds for the troposphere need not be true 
for the stratosphere. In general, the use of an equivalent C0 2 concentration, such as used for 
estimating tropospheric climate change, overestimates the actual stratospheric cooling by the 
non-C0 2 well-mixed greenhouse gases. For realistic mixtures of greenhouse gases, the 
stratospheric cooling is found to be around half of that found using equivalent C0 2 (Wang et al. 
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1991; Shine, 1993). 


The mechanism by which temperatures changes in the stratosphere as a result of constituent 
changes is two-fold. First, the change in the constituent leads to a change in the radiative heating 
field, even if all other conditions are kept fixed. Second, in considering the response of the 
climate system, the timescales for the adjustment of stratospheric temperatures is less than about 
100 days which is much faster than the decadal (or multi-decadal) timescale response of surface 
and lower tropospheric temperatures. This ‘fast’ process of stratospheric temperature change is 
then modified as the troposphere comes into a new equilibrium to the forcing. Thus, for 
example, an increase in C0 2 leads to an initial cooling of the stratosphere. As the surface and 
troposphere warm up due to C0 2 increase, the increased emission from the surface-troposphere 
system contributes a wanning tendency in the lower stratosphere. Forster et al. (1997) have 
shown that this process leads to a significant reduction of the initial cooling in the very lower 
stratosphere for C0 2 ; in the case of CFCs, it leads to a significant heating. When increases in all 
the well-mixed gases over the last century are considered in a ID RCM, the overall result for the 
lower stratosphere is one of an initial radiative cooling, arising due to increase in tropospheric 
opacity which reduces the upwelling longwave flux reaching the stratosphere and its absorption 
there; there is also an increased emission from the stratosphere. The equilibrium result is a 
decrease in lower stratospheric temperatures, with an increased flux divergence of the 
downward longwave beam in the now colder lower stratosphere balancing the increase in the 
flux convergence of the upward beam due to tropospheric warming (Ramaswamy and Bowen, 
1994). Thus, the actual temperature change in the lower stratosphere depends on the degree to 
which the surface-troposphere system equilibrates to the changes in the radiative constituents. 

GCMs have been used to determine the changes in stratospheric climate due to changes in 
the well-mixed greenhouse gases. Fels et al. (1980) obtained a cooling of -10K at 50 km. for 
a doubling of C0 2 . This study also found that FDH results for C0 2 doubling agreed reasonably 
well with the GCM simulations. Subsequently, other GCM experiments of a similar nature 
(e.g., Wang et al., 1991 , Rind et al., 1990) also obtain a substantial cooling of the stratosphere. 
However, the degree of cooling, especially in the polar stratospheres, is model-dependent be- 
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cause of differences m the manner in which physical processes are parameterized in the models 
(Mahlman, 1992; Shindell et al., 1998). GCMs that are used to determine changes in surface 
and tropospheric climates also simulate the stratospheric changes (e.g., those discussed in 
chapter 8 of IPCC, 1996). However, several such models have a coarse vertical resolution in 
the stratosphere which inhibits performing a reliable trends assessment using their results at 
those altitudes. Nonetheless, all models predict a general cooling of the stratosphere due to C0 2 
increases. WMO (1988) estimated a cooling of 0.2-0.3K at 2mb and 0.15K at 10 mb between 
1979/80 and 1985/86 due to C0 2 increases alone. 

A number of ID radiative-photochemical model predictions of the cooling due to increased 
concentrations of greenhouse gases (C0 2 , N 2 0, CH 4 and the CFCs) were reviewed in WMO 
(1988). These models included the temperature feedbacks due to the increased trace gas 
concentrations upon ozone concentration. Our understanding of chemistry at that time was such 
that the important feedbacks were deemed to occur in the upper stratosphere. Briefly, the 
features of these calculations were that peak coolings were estimated to occur between 40 and 
50 km. Simulations for the period since 1960 indicated that the mean cooling rate was about 1.5 
YJ decade, for the period since 1940 it was about 1 K / decade, and for the period from 1850 it 
was about 0.4 YJ decade. There was more spread in model results at 24 km., ranging from 
almost zero to -0.3 YJ decade for the period since 1960. 

Both ID RCMs (e.g., Ramanathan et al., 1985) and 3D GCMs (e.g., Wang et al., 1991) 
indicate that the presence of the well-mixed non-C0 2 trace gases can be expected to yield a 

slight warming of the lower stratosphere. Ramaswamy et al (1996) have compared the 50-100 
hPa temperature change for different time periods for C0 2 alone, and for a realistic mixture of 
gases (Figure 22) using a ID RCM. For the period 1765-1990, the temperature change due to all 
gases is less than 50% of the effect due to C0 2 alone. For more recent periods e.g., since - I960, 
when the relative importance of other gases, and in particular that of the CFCs, has grown 
(Ramanathan et al., 1985), the temperature change due to all well-mixed greenhouse gases is not 
only less than that due to C0 2 alone, but becomes quite small in absolute value. The overall 
cooling effect in the lower stratosphere due to increases in the well-mixed greenhouse gases is to 
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be contrasted with their warming effect on the surface. 


The vertical profile of temperature change in the stratosphere due to recent changes in the 
well-mixed gases is discussed in WMO (1992; see Figure 7-7 of that report). There is a general 
increase of the cooling with height. This is further exemplified by Forster and Shine (1997) who 
show that the FDH-computed increases in the well-mixed greenhouse gases over the decade of 
the 1980s would yield an increasing cooling with height at 40N (-0.7 YJ decade at 35 km.; 
Figure 23). The FDH model of Ramaswamy et al. (1992) considered the well-mixed greenhouse 
gas increases over the 1979-1990 period. The model calculations yield an annual-mean cooling 
that increases with height between -20-50 km. and is fairly uniform throughout the globe 
(Figure 24), with a peak cooling of - 1 YJ decade occurring in the tropics at -50 km. This figure 
also illustrates that the combined well-mixed greenhouse gases’ effect yields a slight but 
characteristic warming near the tropopause, notably in the tropical regions. 


3.3. Stratospheric ozone 

3.3.1 LOWER STRATOSPHERE 

Over the past decade or so, much of the work on stratospheric temperature trends has 
focussed on the effects due to lower stratosphenc ozone change, following the realization of the 
large ozone loss trends in this region, as well as the accumulating evidence for substantial 
temperature changes in this region (WMO: 1990, 1992, 1995). In particular, the large Antarctic 
springtime ozone losses were the first to be examined for their potential temperature effects. 
Subsequently, investigations have been extended to examine the changes initiated by ozone 
depletion m the global lower stratosphere. 

As a general demonstration of the sensitivity of the global lower stratosphere to changes in 
ozone. Figure 25 shows the change in lower stratospheric and tropospheric temperatures when 
the entire ozone in the 70-250 hPa region is removed. Temperature decreases of up to about 4K 
are obtained. Note that the region in the vicinity of the tropopause is most substantially affected, 
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consistent with the known radiative sensitivity and large radiative damping time of this region 
(Fels, 1982; Kiehl and Solomon, 1986). A warming is seen above the cooling which is in part 
due to more upwelling thermal infrared radiation reaching this region in the absence of ozone in 
the 70-250 hPa layer, and in part due to dynamical changes; these will be elaborated upon later. 

With regards to actual decadal trends, Miller et al. (1992) estimated a cooling of the lower 
stratosphere by about 0.3K/ decade over the 1970-1986 period based on rawinsonde 
observations. By performing a calculation of temperature changes using observed ozone trends, 
they inferred a substantial role due to ozone losses in the observed temperature trend. Since 
ozone changes vary greatly with latitude, there is a need to proceed beyond global-means and 
perform model calculations which resolve the latitudinal variations. Locally, the stratosphere is 
not in radiative equilibrium, and the temperature is determined by both radiative and dynamical 
processes. However, by considering the stratosphere to be in a local radiative-dynamical 
equilibrium, with the dynamical heating rates fixed over some specific time scale (e.g., season), 
the FDH concept in effect has been extended to derive the local columnar temperature changes. 
Shine ( 1 986) demonstrated that a large cooling would occur in the Antarctic winter/ spring 
lower stratosphere owing to the ‘ozone hole’. 

McCormack and Hood ( 1 994) presented the first attempt to match the observed seasonal and 
latitudinal variation of lower stratospheric temperature with that predicted using an FDH model 
using imposed ozone changes based on observations. The resemblance between the size and 
pattern of the model temperature changes, with peak changes in the winter/ spring of northern 
mid-latitudes and spring of southern high latitudes, and MSU satellite observations (Randel and 
Cobb, 1994) was found to be encouraging. As discussed below, this work has been taken further 
by using GCMs, which strengthens the model-observation comparisons. 

One particular item of interest in the case of FDH calculations (e.g., McCormack and Hood, 
1994) is that the peak cooling in the Antarctic springtime tends to occur about one month earlier 
than in the observations. This occurs because the standard FDH approximation does not account 
for the time period during which the perturbation in the constituent (in this case ozone) persists. 
In the Antarctic spring time, typical radiative timescales are in excess of a month and yet the 
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timescale of the most marked ozone change is only one month. Hence the atmosphere does not 
have time to fully equilibrate to the ozone-induced radiative perturbation (Shine 1986); the 
standard FDH approximation does not account for this, and assumes that the temperature 
adjustment occurs essentially instantaneously. Figure 26 from Forster et al. (1997) shows the 
temperature change at 80S and 18.5 km. using seasonally-varying ozone changes m the context 
of the standard FDH approach, and a modified FDH approach (termed the seasonally-evolving 
FDH or SEFDH) which accounts for the timescale of the perturbation. The SEFDH result is 
obtained by time-marching the radiation calculations through the seasons. As can be seen, the 
maximum cooling in the FDH occurs earlier in the austral spring and is significantly larger than 
the SEFDH-derived maximum cooling Also notable is that the FDH model does not retain a 
‘memory’ of perturbations at other times of year. An ozone loss in October only would have no 
effect on temperatures at other times of year. The SEFDH model does retain a memory of the 
radiative perturbations, and the increased cooling which persists from late spring to late summer 
in the SEFDH calculations is because the model temperatures have not yet recovered from the 
large loss of ozone in the springtime. 

Early GCM works investigating the effects of ozone depletion did so assuming uniform 
stratospheric ozone losses, with a view towards diagnosing in a simple manner the radiative- 
dynamical response of the stratosphere (e.g., Fels et al., 1980; Kiehl and Boville, 1988). Al- 
though the depletion profiles employed did not quite represent those observed subsequently in 
the lower stratosphere, these studies have shed insights on the mechanism of the stratospheric 
response to ozone changes there Recently, Christiansen et al. ( 1 997) have further examined the 
nature of the stratospheric response to different types of idealized ozone perturbations. 

With the steady build-up of knowledge since the early-to-mid 1980s on the spatial and tem- 
poral distribution of ozone losses in the lower stratosphere, GCM experiments of two types 
have been carried out to determine the resulting temperature changes in the stratosphere. One 
type of investigation concerns attempts to simulate the interactive radiative-chemistry -dynam- 
ical changes due to ozone losses in the Antarctic polar region during springtime e.g., Cariolle 
et al. (1990) and Prather et al. (1990) found a substantial lower stratospheric cooling during the 
duration of the ozone hole. Mahlman et al. (1994) found a statistically significant cooling in the 
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Antarctic spring due to the ozone loss. They also found a significant warming of the altitude 
region above the altitudes of ozone loss, consistent with the observations (Randel, 1988; Figure 
15, 16). Recent studies have extended the interactive simulations to the Arctic region as well. 
Austin et al. (1992) concluded from their investigation that the stratosphere can be expected to 
be cooler owing to both long-term C0 2 increases and ozone decreases, which raises the possi- 
bility in the future of Arctic ‘ozone holes’ and a strong positive feedback effect involving ozone 
depletion and lower stratospheric cooling. Shindell et al. (1998) have also examined the effects 
of increased well-mixed greenhouse gases and ozone changes upon the global stratospheric 
temperatures over the next few decades (see Chapter 12). Their results reiterate the sensitivity 
of the polar springtime temperatures to greenhouse gas increases and the feedback effect 
caused by ozone losses, with the actual magnitudes subject to details of model representations 
of various physical processes. 

Another type of GCM investigation has consisted of imposing the observed ozone losses and 
then determining the GCM response, without any considerations of the chemical and dynamical 
processes affecting ozone distributions. Again, early studies began with investigations of the Ant- 
arctic springtime lower stratosphere owing to the large observed losses there. A substantial lead-in 
into the 3D aspects was provided by the 2D study of Chipperfield and Pyle (1988). They found a 
large cooling (>9K) at -70 hPa during October, just like the earlier FDH result (Shine, 1986). In 
addition, they found a warming (up to -6K) in the upper stratosphere. The GCM simulation of Kie- 
hl et al. (1988) yielded a large cooling (~5K) by October end, which further substantiated the FDH 
inference. This GCM study also found a warming in the upper stratosphere (up to ~4K at 10 hPa). 

Employing the observed global lower stratospheric ozone depletion, Hansen et al. (1995) dem- 
onstrate a good agreement in the magnitude of the GCM-simulated cooling of the global, annual- 
mean lower stratosphere with the observed MSU satellite and radiosonde temperature trends (Fig. 
27). This study along with others also indicates that the cooling due to ozone loss overwhelms the 
temperature change resulting from changes in concentration of the well-mixed greenhouse gases 
in the 50-100 hPa (-16-21 km.) region (see Figures 22, 23 and 27). 

The differences between FDH and GCM simulations of the effects of ozone losses have been 
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investigated by Fels et al. (1980) and Kiehl and Boville (1988) using idealized stratospheric ozone 
perturbations. For annual-mean model simulations, Fels et al. (1980) concluded that the FDH re- 
sponse to uniform percentage reductions tended to be comparable to the GCM solutions, except at 
the tropical tropopause and tropical mesosphere. For “perpetual January” simulations, Kiehl and 
Boville ( 1 988) found that the uniform reduction experiments were not well represented by the FDH 
calculations. However, more realistic ozone scenarios in which the losses increase from equator to 
the poles yielded FDH results comparable to the GCM. This point is substantiated by the “perpet- 
ual January” simulations of Christiansen et al. (1997), who performed uniform ozone reduction and 
large lower stratospheric depletion experiments. The main dynamical effect in their perturbation 
experiments was a weakening of the diabatic meridional circulation accompanied by a latitudinal 
smoothing of the temperature response. 

Using the satellite-observed global lower stratospheric ozone losses over the -1979-1991 peri- 
od (i.e., just prior to the Pinatubo volcanic eruption), a comparison of the resulting radiative and 
radiative-dynamical solutions for the stratospheric temperature changes can be obtained. This is il- 
lustrated here by comparing the results of the FDH and GCM simulations peformed, respectively, 
by Ramaswamy et al. (1992) and (1996) using an identical model framework (Figure 28). This 
comparison enables a delineation of the role of dynamical influences on the temperature changes 
caused by the observed ozone depletion. The GCM result, like the FDH, indicates a cooling of the 
lower stratosphere, but there are distinct differences due to dynamical changes. In the mid-to-high 
southern latitudes, there is less cooling in the GCM. In the northern hemisphere, the midlatitudes 
are less cold in the GCM, but the high latitudes are more so relative to the FDH result, again a con- 
sequence of the dynamical changes in the model. In the GCM, there is a cooling even in those re- 
gions where there are no ozone losses imposed e.g., the lower stratosphere equatorward of 15 
degrees. A warming occurs above the region of cooling, particularly noticeable in the southern 
hemisphere, similar to the results obtained by other GCM studies (Kiehl et al., 1988; Mahlman et 
al., 1994; Shindell et al., 1998). The dynamical changes (see also Mahlman et al., 1994) consist of 
an induced net rising motion in the tropics and a compressional heating of the middle stratosphere 
at the higher latitudes. The annual-mean response is statistically significant between -13 and 21 
km. in the -20 to 50 degree latitude belt (Ramaswamy et al., 1996). The changes at high latitudes 
(>60 degrees) fail the significance test because of large interannual variability in those regions. The 
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warming above the lower stratospheric regions in both hemispheres is reasonably similar to obser- 
vations (Randel, 1988; Figures 15, 16). On the basis of the GCMs-observations comparisons, a 
principal factor for the observed warming could be the radiative-dynamical feedbacks involving 
ozone depletion in the lower stratosphere. 

Hansen et al. (1993) show that the zonal-mean patterns of GCM-simulated lower stratospheric 
temperature change due to imposed ozone losses correspond well with observed changes. In a 
study analogous to the FDH model -observation comparison discussed earlier, Ramaswamy et al. 
(1996) have compared the latitude-month trend pattern of the decadal (period: -1979-1990) tem- 
perature change and its statistical significance, as simulated by a GCM (Figure 29a) in the altitude 
region of the observed lower stratospheric ozone change (tropopause to -7 km. above), with that 
derived from satellite observations (Figure 29b) of the lower stratosphere for the same period (Ran- 
del and Cobb, 1994). The match of the latitudinal-month pattern of cooling with observations bears 
a fair resemblance to the FDH-based study of McCormack and Hood (1994). In the midlatitudes, 
both panels illustrate a cooling from -January to October in the northern hemisphere and from 
-September to July in the southern hemisphere. The cooling in the midlatitudes of the northern 
hemisphere from -December to July, and in the southern hemisphere from -December to May, are 
statistically significant in both model and observation. Comparing with Figure 13, it is apparent 
that the observed space-time domain of statistically significant cooling is also dependent on the end 
year chosen for the analysis. Near the poles, both the simulation and observation exhibit relatively 
large magnitude of cooling during winter and spring. The simulated cooling in the Antarctic is 
highly significant during the austral spring (period of the ‘ozone hole’), consistent with observa- 
tion. The springtime cooling in the Arctic does not show a high significance owing to a large dy- 
namical variability there. The simulated cooling in the tropics, which arises as a result of changes 
in circulation and is absent in FDH, is not significant for most of the year owing to small tempera- 
ture changes there. There exists quantitative differences between the simulated and observed 
trends. In addition, there is less variability in the model compared to observations. 

In general, the observed ozone losses introduce a nonuniform space-time cooling. This is in 
contrast to well-mixed greenhouse gas effects which tend to yield a more uniform lower strato- 
spheric cooling with season (not shown). Thus, the observed seasonal trend (Figure 29b), which 


44 



exhibits a spatial dependence, is unlikely to be due to solely increases in well-mixed greenhouse 
gases. Note that, in the context of the GCM results, it is not possible to associate specific seasonal 
cooling at any location with the corresponding seasonal ozone loss. 

Some models that have been principally used to study tropospheric climate changes due to 
greenhouse gases (e.g., IPCC, 1996) have also considered the lower stratospheric ozone deple- 
tion. However, as noted earlier, several of these climate models have a poor vertical resolution 
in the stratosphere which affects the accuracy of the calculated temperature changes. Neverthe- 
less, the cooling predicted by these models using observed ozone loss in the lower stratosphere 
is generally consistent with the observed global-mean cooling trend, typically 0.5 K/ decade. 

In addition to lower stratospheric ozone, model calculations illustrate that tropospheric 
ozone change could also have a smaller but non-negligible effect on stratospheric temperatures 
(Ramaswamy and Bowen, 1994). Since pre-industrial times, tropospheric ozone may have in- 
creased by -50% or more in some regions in the Northern Hemisphere (IPCC 1995; Bemtsen 
et al., 1997). As a result, this may have contributed to a substantial cooling tendency in the low- 
er stratosphere in specific regions. 

In a diagnostic study, Fortuin and Kelder (1996) discuss relationships between ozone and 
temperature changes by analyzing concurrent measurements of these two quantities made at 
eight stations over the past two decades. They find that both vertical displacement and radiative 
adjustment to trace gas concentrations need to be considered in order to explain the observed 
ozone and temperature variations at the various locations. 

3.3.2 SENSITIVITIES RELATED TO OZONE CHANGE 

One of the most important factors affecting the temperature change due to ozone is the 
precise vertical profile of the ozone changes in the stratosphere, especially the region near the 
tropopause, a problem that has been reiterated for a number of years now (WMO: 1986, 1992, 
1995). Additionally, there are some uncertainties about the middle and upper stratospheric 
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ozone changes that, too, impact in a non-negligible manner upon the temperature changes in 
those regions. 

The sensitivity of the temperature change to the details of the vertical profile of ozone 
change has been discussed by Forster and Shine (1997). Figure 30 shows the latitude-height 
profile of the FDH-derived temperature change using a simplified vertical profile of ozone 
change (where a constant percentage of ozone is removed within 7 km. of the tropopause, with 
the total column ozone change derived from the SBUV trends reported in WMO, 1995), and one 
where the vertical profile of trends derived from SAGE data is used. The calculations are for the 
period since 1979. The vertical profile of the cooling obtained due to the two ozone profiles at 
40N is compared with other species in Figure 23. Generally speaking, the cooling in the lower 
stratosphere due to either profile is strong and outweighs contributions estimated due to other 
possible species’ changes. Although the general vertical profile of the ozone-induced coolings 
remains the same, peaking at ~ 1 YJ decade in some regions, the latitudinal variation is quite 
different. The SBUV trends give only small ozone changes in low latitudes, and the cooling is 
concentrated in the extratropical lower stratosphere. The SAGE trends show large changes in the 
tropical lower stratosphere (despite modest total column changes) and these result in the peak 
cooling in the tropics. It is noted that the SPARC Ozone Trends Panel has recently assessed the 
status of knowledge about the vertical trends of ozone changes. According to that panel and 
Cunnold et al. (1996a,b), the large trends of SAGE I/D in the 15-20 km. are unrealistic, with 
these data being most reliable above 20 km. 

One further important aspect of Figure 30 concerns the trends in the middle stratosphere. 
Using the simplified ozone change profile, there is only heating at -30 km., most marked in 
southern high-latitudes where the increase in the upwelling thermal infrared radiation, as a result 
of the removal of ozone in the lower stratosphere, leads to an increase in the absorbed radiation. 
Using the SAGE vertical profiles, which are probably more reliable at these heights, there are 
coolings at almost all latitudes, reaching 0.4 K / decade. Hence, whilst the ozone loss is probably 
most important in the lower stratosphere, the losses at these stratospheric altitudes (-30 km.) 
remain substantial, and could lead to trends that are of the same magnitude approximately as the 
temperature trends, over the same period, due to increases in well-mixed greenhouse gases (see 


46 



Figure 24). 


Hansen et al. (1997) have examined the sensitivity of the lower stratospheric temperature 
changes due to ozone depletion by conducting various GCM experiments with varying profiles 
and amounts of ozone losses (their result for the specific case when depletion of ozone occurs 
over the 70-250 hPa region is illustrated by Figure 25). They find that the magnitude of tem- 
perature change near the tropopause is strongly dependent on the vertical profile of the ozone 
loss. The degree to which the temperature change is affected is also dependent on the latitude 
under consideration. 

The GCM sensitivity experiments by Hansen et al. (1997) mentioned above also demon- 
strate that the observed ozone losses not only cool the lower stratosphere but also affect the ver- 
tical profile of temperature change below the region of ozone losses. The cross-over of 
temperature change from a warming in the troposphere to a cooling in the stratosphere occurs 
at a lower altitude when stratospheric ozone depletion is considered in the models, in contrast 
to when only the changes in the well-mixed gases are considered. Santer et al. (1996) have 
shown how the inclusion of the ozone loss in the lower stratosphere not only alters the modeled 
vertical profile of the atmospheric temperature change due to anthropogenic species, but is a 
critical element in the attribution of the observed temperature profile change to human activity. 
Tett et al. ( 1 996) and Folland et al. ( 1 998) further point out that the inclusion of the ower strato- 
spheric ozone depletion is an important anthropogenic factor for the vertical profile of temper- 
ature changes near the tropopause, and thereby could affect the lower altitudes and surface, too. 

Ozone changes also affect higher regions of the stratosphere. Fels et al. (1980) and Rind et 
al. (1990, 1998) point out that, while the stratosphere generally cools in response to the carbon 
dioxide increases, the cooling in the middle and upper stratosphere is reduced (by -10% at 50 
km.) when ozone is allowed to respond photochemically. FDH estimates using the SAGE ver- 
tical profile of temperature change indicate that the positive (or negative) temperature change 
generally follows the positive (or negative) changes in ozone. Thus, the result from the calcu- 
lations performed by Schwarzkopf and Ramaswamy (1993) suggests that ozone losses over the 
1980-1990 period have caused a cooling of -0.3 K/decade at -40 km. This is to be contrasted 
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with the well-mixed greenhouse gas-induced cooling of about 0.8 K/ decade at these altitudes 
(Figures 23, 24). 

As mentioned in section 3.3.1, Shindell et al. (1998) simulate the changes in the strato- 
sphere due to secular trends in greenhouse gases and ozone changes using an interactive radi- 
ative-dynamical-chemical GCM (see also Chapter 12). Their results indicate that the cooling 
of the entire stratosphere due to the trace gas changes leads to a significant expansion of the 
area of ozone loss in the polar regions. Besides the expansion of the Antarctic polar ozone loss 
area, there also occurs a pronounced depletion of the Arctic springtime ozone. The rather dra- 
matic modeled cooling is attributed to changes in planetary wave characteristics. In general, the 
quantitative nature of such changes in GCMs depend crucially on the model parameterizations. 

3.4 Aerosols 

As stated in section 2.6. 1 , major volcanic eruptions result in temporary increases of strato- 
spheric aerosol concentrations which, in turn, cause transient warmings of the lower strato- 
sphere (WMO: 1988, 1992). One-dimensional radiative-convective, two-dimensional and 
three-dimensional models have all been used to estimate the effects on stratospheric tempera- 
tures (e.g., Hansen et al., 1978; Harshvardhan, 1979; Pitan et al., 1987; Hansen et al., 1997b). 
An important research outcome in recent years has been the fact that simulations of the climatic 
effect of the Pinatubo eruption (Hansen et al., 1993; Hansen et al., 1997b) reproduce the ob- 
served lower stratospheric warming (refer to section 2.6.1; Christy, 1995; Randel et al., 1995; 
Angell, 1997; also Figures 2 (bottom panel) and 3) - the peak and its duration - remarkably well 
(IPCC, 1995). Although temporary, the volcanic aerosol-induced lower stratospheric warming 
is in the opposite sense of the cooling due to ozone depletion, and the cooling due to the com- 
bined well-mixed greenhouse gas increases (Figures 23, 24, 27). 

Rind et al. (1992) investigated the role of volcanic effects on the stratosphere in GCM ex- 
periments with and without sea-surface temperature changes. Their results show that, in re- 
sponse to the volcanic injections, the lower stratospheric warmed, while the upper stratosphere 
cooled and the mesosphere warmed at mid-and-high latitudes. These responses are similar 
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qualitatively to the OHP lidar observations. The occurrence of mesospheric warming is indi- 
rectly supported by a decrease in the occurrence of noctilucent clouds in the summer of 1992 
(Zalcik, 1993). 

Several 2D radiative-dynamical-photochemical models (Brasseur and Granier, 1992; Tie et 
al., 1994, Rosenfield et al., 1997), with different degrees of complexity in the representation of 
physical and chemical processes, have been deployed to study the temperature, circulation and 
ozone changes due to the increase in stratospheric aerosol concentrations following the Mt. Pi- 
natubo volcanic eruption. Additionally, Eluszkiewicz et al. (1997) have used the UARS mea- 
surements to investigate the stratospheric heating rate sensitivity and the resulting residual 
circulation. These studies indicate a warming of the tropical lower stratosphere that is approx- 
imately consistent with observations. This heating is found to be sensitive to tropospheric 
clouds’ presence and to the aerosol vertical extinction profiles which corroborates the sensitiv- 
ity results of WMO (1988). The perturbed heating in the tropics results in an upwelling and an 
altered circulation, a feature seen in model simulations and diagnosed from observations. Thus, 
the transient volcanic aerosol-induced warming has the potential to alter atmospheric dynamics 
and could thereby affect the surface-troposphere system (Hansen et al., 1997b). Ozone deple- 
tions are obtained by the 2D model as a result of the heterogeneous processes involving aero- 
sols, changes in photolysis rates and changes in the meridional circulation (Tie et al., 1994). 
The ozone loss at the high latitudes yields a cooling there (Rosenfield et al., 1997) which could 
have contributed to the prolonged duration of the colder-than-normal temperatures in the Arctic 
region (section 2.4; Figure 4, 14(top)). 

Although both El Chichon and Pinatubo eruptions have resulted in a large transient 
warming, these are estimated to have only a small direct effect on the decadal temperature 
trends in the lower stratosphere, with even smaller effects higher up. However, an indirect effect 
through aerosol-induced, chemically-catalyzed ozone loss, and the subsequent ozone radiative 
feedback, can yield an enhanced global cooling trend (Solomon et al., 1996). Figures 2 (bottom 
panel) and 3 illustrate the colder global-mean lower stratospheric temperatures now than in the 
pre-Pinatubo period (see also Chapter 7). Thus, the volcanic events could have had an important 
bearing on the lower stratospheric temperature trend observed over the decade of the 1980s and 
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early 1990s. 


Ramaswamy and Bowen (1994) show that changes in tropospheric aerosols can affect lower 
stratospheric temperatures. While the tropospheric aerosols assumed in that study acted to 
oppose the greenhouse warming in the troposphere, they enhanced the cooling in the 
stratosphere by reducing the upwelhng thermal infrared radiation coming from a cooler 
troposphere. It should be noted that there is considerable uncertainty in the distribution and 
optical properties of aerosols; in addition, they may have a substantial effect on the optical 
properties of clouds (see IPCC: 1995, 1996) which would also affect the upwelling longwave 
flux from the troposphere. Thus, it is not possible as yet to be quantitative about the present-day 
tropospheric aerosols’ effect on stratospheric temperatures. 

3.5 Water vapor 

Stratospheric water vapor is an important radiatively-active constituent, and it is expected to 
increase in concentration due to increased methane oxidation. It could also be influenced by 
changes in tropospheric water vapor concentration and changes in the transport of water vapor 
across the tropopause. Oltmans and Hofmann (1995) reported significant trends of 0.5%/year in 
water vapor over Boulder Colorado from balloonsonde observations between 1981 and 1994. It 
is not clear whether these changes are occurring at other locations and, more importantly, on a 
global scale. Forster and Shine (1997) have used a FDH model to calculate the possible impact 
of these changes. At 40N and between 15 and 20 km., the cooling due to water vapor changes 
exceeds 0.2 K / decade (Figure 23). This is much larger than the effect of changes in the well- 
mixed greenhouse gases at the same altitude (typically 0.1 K/ decade; Figures 23, 24); 
depending on the vertical profile of ozone loss adopted, the water vapor change causes a cooling 
of around 20-30% of the effect of the ozone change. Recent observational analyses (Nedoluha et 
al., 1998) suggest an increase in H 2 0 of 1-2% or more per year (between 1991 and 1996) in the 
middle and upper stratosphere and mesosphere. If sustained for a decade, the heating rate 
change would be about 0.05 YJ day. Combining this with a radiative damping time of ~10 days 
(Fels, 1982; Kiehl and Solomon, 1986) would yield a cooling of 0.5 YJ decade at 50 km. This 
value would be about half of that for the well-mixed greenhouse gases (section 3.2; Figure 24), 
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and could exceed that due to ozone change at that altitude (section 3.3.2). 


3.6 Other (solar cycle, QBO) 

All the above changes in constituent changes have been ones that have taken place in the 
stratosphere itself. External changes, such as changes in the incoming solar radiation, reflected 
solar radiation from the troposphere, and in the emitted thermal infrared radiation can also 
impact stratospheric temperatures. 

Changes in the output of the Sun can clearly influence the amount of UV and visible 
radiation absorbed by the stratosphere. In the absence of changes in ozone, WMO (1988) 
radiative-convective model calculations showed that, over the course of a solar cycle the 
temperature at the stratopause can vary by about IK, but by 20 km. the changes are less than 
0.1K. If planetar}' wave propagation plays a role in the way the solar cycle influences the 
stratosphere as observations appear to indicate, then model simulations need to adequately 
account for these dynamical processes. Several numerical model simulations (Garcia et al,1984, 
Brasseur, 1993, Fleming et al. 1995, McCormack and Hood, 1996) indicate that the maximum 
of solar effect should be situated at the stratopause at all latitudes, with a value of around 1.5- 
2K / cycle in the tropics. Figure 3 1 illustrates results from two models for the equatorial latitudes 
and compares them with rocket and satellite (see also Figure 19) observations. The models yield 
amplitudes of ~2K at the tropical stratopause. The modeled temperature change is roughly 
consistent with observations, but substantial differences remain; the same models are less 
successful when compared with observations at middle and high latitudes. In fact, any potential 
dynamical effects associated with the observed solar cycle signature (section 2.6.2) are likely 
not being captured adequately by interactive radiative-chemical-dynamical models as yet. 

Haigh (1996) has simulated the possible effect of solar-induced stratospheric ozone and 
temperature variations on the tropospheric Hadley circulation and jet streams using a GCM. 
The results suggest that changes in stratospheric ozone may provide a mechanism whereby 
small changes in solar irradiance can cause a measurable impact on climate. Haigh (1996) pro- 
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posed that changes in the thermal structure of the lower stratosphere could cause a strengthen- 
ing of the low latitude upper tropospheric easterly winds. Haigh (1998) further suggests that 
this easterly acceleration can come about els a result of a lowering of the tropopause and decel- 
eration of the poleward meridional winds in the upper portions of the Hadley cells. More ob- 
servational diagnostics work and model investigations are needed to confirm these hypotheses. 

Balachandran and Rind (1995) studied the influence of UV variability and QBO using a 
GCM that emphasizes the middle atmosphere. They found that changes in UV tended to alter 
the radiative heating with an accompanying influence on the winter polar jet, though it must be 
noted that the large variability in this quantity inhibits a statistically significant and robust de- 
termination. The high latitude response to UV is affected by alteration in planetEiry wave activ- 
ity, with the QBO exerting a modulating effect. The influence of solar change on the 
temperature yielded changes of opposite sign at the stratopause between latitudes above and 
below 60N. Such changes are quite similar to those seen in the overall stratosphere from the 
SSU satellite analyses (Figure 19). 


4. CHANGES IN TRACE SPECIES AND OBSERVED TEMPERATURE TRENDS 

4.1 Lower stratosphere 

As far as the global lower stratosphere is concerned, there is now a firm documentation of 
the changes in ozone and greenhouse gases. The importance of ozone depletion relative to that 
due to changes in other greenhouse gases that are well-mixed (C0 2 , CH 4 , N 2 0, CFCs) have 
been evaluated by several studies with various types of models (ID to 3D). Miller et al. (1992) 
and Hansen et al. (1995) demonstrate that the global-mean lower stratospheric temperature 
change in the 1980s can be explained only when ozone changes are considered (Figure 27). Ra- 
maswamy et al. (1996) show that the global, annual-mean GCM temperature change due to the 
decadal ozone losses in the -50-100 hPa (-16-21 km.) lower stratospheric region is much 
greater than that due to increases in C0 2 only and all well-mixed greenhouse gases taken to- 
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gether (Figure 22). It is thus encouraging that all model studies, ranging from FDH and RCMs 
to GCMs yield similar findings (section 3 3.1). The global-mean decadal cooling in the 1980s 
due to ozone is estimated to be -0.5 to 0.6K which is comparable to the reported decadal trends 
from observations (section 2.3.2; Figure 27). The well-mixed gases’ effect is less than one- 
fourth that due to ozone depletion. While the increase in C0 2 alone since 1765 yields a cooling 
of -0.3K, inclusion of the other well-mixed gases, which together tend to warm the tropopause 
region, yields -0.15K. It is thus clear that the computed 1979 to 1990 ozone effect on lower 
stratospheric temperature outweighs the effects of changes in other gases not only over the last 
10 to 30 years, but also over the past two centuries. In the global, annual-mean, it can be strong- 
ly declared that the observed ozone depletion in the lower stratosphere is the dominant radiative 
cause of the observed global-mean cooling. [Chapter 7, WMO (1998) presents a theoretical ba- 
sis for how a long-term global, annual-mean trend in the stratosphere must necessarily be as- 
sociated with radiative perturbations ] 

As discussed in section 3.3.2, uncertainties arise in the simulation owing to incomplete ob- 
servational knowledge of the vertical profile of global ozone loss near the tropopause, including 
that in the tropical areas. While more thorough altitudinal measurements of ozone loss would 
lead to more precise simulation of temperature change, with cooling extending to perhaps even 
higher altitudes (e.g., springtime southern polar latitudes), the lower stratosphere region, taken 
as a whole, can be expected to cool notably given the magnitude of the ozone losses observed. 
As a principal conclusion from all investigations thus far, whether ozone changes are pre- 
scribed or determined self-consistently within a model, the global lower stratospheric region, 
especially the mid-to-high latitudes, cools in a significant manner in simulations for the de- 
cades of the 1980s and 1990s 

In addition to the above, a principal feature from especially the GCM studies is the reason- 
ably good correspondence of the zonal-mean lower stratospheric cooling trends since -1979 
with satellite and radiosonde records (e g., Hansen et al., 1993). It is concluded that the reason- 
able consistency of the simulated cooling pattern and magnitudes with that observed, including 
the regimes of statistically significant changes, coupled with the high correlations noted be- 
tween observed temperature changes and ozone losses (section 3.3.1), confirm the notion that 
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ozone depletion has caused a substantial spatially-and-seasonally-dependent effect in the lower 
stratosphere over the past decade. For example. Figure 29 highlights the model-observation 
consistency with regard to magnitude and statistical significance in the midlatitudes of both 
hemispheres during the first half of the year, and during the Antarctic springtime. Although the 
attribution of the observed temperature trends to the observed ozone depletion in the lower 
stratosphere is strong in the global, annual-mean sense, the spatial and temporal aspects de- 
mand more circumspection. However, no other cause besides ozone depletion has been shown 
as yet to yield a latitude-month fingerprint such as that seen in the observations (Figure 29b). 
Thus, in the zonal-mean, seasonal sense, it can be stated that ozone is identified as an important 
causal factor of lower stratospheric temperature change. 

Possible secular changes in other radiatively-active species are estimated to contribute 
smaller decadal effects than the stratospheric ozone loss. Information on decadal changes in 
global water vapor (section 3 5) and clouds is insufficient to estimate their influence precisely; 
there is no information at present to suggest that their effects could be as dominant as that due 
to the stratospheric ozone loss. Although volcanic aerosols can have substantial impact over the 
1-2 years that they are present in the lower stratosphere, their effect on the past decade’s tem- 
perature trend has probably been small compared to ozone (sections 2.6.1, 3.4). There is little 
evidence to suggest that forcings from the troposphere (e.g., sea-surface temperature changes; 
section 2.6.5) or natural climate variability or solar cycle (sections 2.6.2, 3.6) have significantly 
influenced the global lower stratospheric temperature change over the past 1.5 decades, al- 
though, in the absence of rigorous long-term observations, a precise estimate of their 
contributions cannot be obtained. It is noted that some ozone loss has been reported for the 
1970s, too, (Lacis et al., 1990) which would have contributed to the small observed cooling 
during that decade (e.g., Figure 3(top), 4). However, this contribution is not likely to have been 
as much as in the 1980s and 1990s, since the ozone losses for the earlier decades are concluded 
to have been never as high as those in recent ones (Bojkov and Fioletov, 1995). 

There is a scarcity of knowledge on the low frequency variability of the stratosphere and its 
causes - from either observations or models. Unlike surface temperature measurements which 
span multi-decades, those for the stratosphere are available only from about late 1950s, and the 
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continuous record is available only at a few locations in the NH. This makes it difficult to assess 
accurately the low frequency variability. While global coverage has become possible with the 
MSU and SSU satellites since 1979, the time period available to-date is too short to assess any- 
thing beyond mterannual variability; certainly, rigorous decadal-scale variability analysis will 
not be possible until data for a few more decades become available. The model simulations to- 
date suggest an interannual variability in some features that bear some resemblance to that ob- 
served, but there are also features that the existing models either cannot reproduce or fail to 
mimic the observations (e.g., Hamilton et al., 1995). A prominent uncertainty arises due to the 
lack of a proper simulation of the polar wintertime and winter-to-spring transitional tempera- 
tures (including sudden warmings) from first principles. The usual method to reproduce obser- 
vations is to “tune” the model in some manner e.g., gravity-wave drag. The quantitative effect 
that this has on the fidelity of the simulation of trends and variability remains to be determined. 

The cold winters/ springs occurring in the lower stratosphere of the Arctic in recent times 
(section 2.4) have drawn considerable attention, including a search for the possible cause/s. The 
Arctic temperatures can vary substantially in magnitude on interannual time scales (see Figure 
14(top)), thus demanding caution in the inference and attribution of long-term trends there. 
Baldwin et al. (1994) and Perlwitz and Graf (1995) document the existence of a coupling be- 
tween the winter polar vortex and middle troposphere. Thompson and Wallace (1998) show co- 
herence between hemispheric wintertime surface air temperature and pressure and 50 hPa 
temperatures. The studies of Manney et al. (1996) and Coy et al. (1997) indicate that extreme 
low temperatures in spring are associated with record low planetary wave activity and an in- 
tense vortex. The above studies suggest a possible shift in Arctic meteorological conditions in 
early-to-mid 1990s, one that is linked to tropospheric circulation and variability in planetary 
wave forcing. The notion that this could be related in some manner to greenhouse gas-induced 
warming of the lower atmosphere (Chapter 12) remains to be conclusively demonstrated. 

Indeed, the effect of tropospheric climate change due to changes in trace gases and aerosols, 
equator-to-pole tropospheric temperature gradient, waves propagating into the stratosphere, 
and the ensuing effects on the radiative-dynamical-chemical stratospheric equilibrium are not 
well understood in a quantitative manner. For example, different GCMs suggest substantially 
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different manner of changes in the characteristics of the planetary wave activity due to increase 
of C0 2 , which would, in tum, impact the radiative-dynamical interactions and the magnitude 
of stratospheric temperature changes (Fels et al., 1980; Fels, 1985; Rind et al., 1990; Mahlman 
et al., 1992; Graf et al., 1995; Shindell et al., 1998; Rind et al., 1998). 


4.2 Middle and upper stratosphere 

Unlike the case for the lower stratosphere, the trends estimated from the different observa- 
tional platforms for the middle and upper stratosphere are not as robust. The satellite, lidar and 
rocket data, although having a consistency of pattern above ~50 hPa (e.g.. Figures 9 and 12), 
do not exhibit the same degree of coherency that exists with respect to both magnitude and sta- 
tistical significance for the different datasets of the lower stratosphere (Figures 5, 6). 

In the middle and upper stratosphere, model results suggest that the increases in the well- 
mixed greenhouse gases and changes in ozone will contribute to temperature changes. Figure 
24 shows the zonal, annual-mean temperature trend due to the well-mixed greenhouse gas in- 
creases in the 1980s. The overall picture in the annual-mean is one of cooling from the lower 
to the upper stratosphere (see also Figure 23). This cooling in the middle stratosphere due to 
the well-mixed gases can be expected to be enhanced in FDH calculations that also consider 
ozone losses in the middle stratosphere; the latter is estimated to yield about a 0.3 K/decade 
cooling using the SAGE depletion for the 1980s period (section 3.3.2). The computed vertical 
profile bears a qualitative similarity to the observations (e.g., Figures 7, 9) with regard to the 
cooling of the entire stratosphere. This reaffirms the secular cooling trend due to greenhouse 
gas increases inferred for the stratosphere from shorter records (WMO: 1988, 1990). 

However, at altitudes above the lower stratosphere, there are major quantitative differences 
between the modeled and observed cooling. The FDH simulated cooling increases with height 
when only the well-mixed gases are considered (Figure 24), whereas the observations indicate 
a rather uniform trend between 20-35 km. (with perhaps even a slight reduction at -30-35 km.; 
see Figures 7 to 1 2). It must be noted here that the GCM’s simulation of the effects due to ozone 
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loss involve a dynamical change (not evident in FDH) that causes some warming above the lo- 
cation of the cooling in the lower stratosphere (section 3.3.1 and Figure 28). This GCM-simu- 
lated warming, which may not be statistically significant, occurs mainly in winter/ springtime, 
and is qualitatively consistent with observations during that season (Figures 15, 16). Consider- 
ation of this warming could retard somewhat the rapid increase of the well-mixed greenhouse 
gas-induced cooling with height seen in Figure 24 (e.g., above 25 km.) and enable a full radi- 
ative-dynamical (i.e., GCM) computation to become more consistent with observations than 
the FDH results. 

Additionally, the magnitude of the modeled cooling in the upper stratosphere is less than 
that observed e.g., at -45 km., the modeled cooling is about 1 K / decade due to the well-mixed 
greenhouse gases and about 0.3K/ decade due to ozone, while the observed cooling is greater 
than 1.5K/ decade (e g., at 45 km. in Figures 9, 12). Some of this bias could be due to water 
vapor whose decadal trend in the 1980s globally is not known. Recent satellite data suggests 
an upward trend over the past five years which would add to the cooling trend computed for the 
upper stratosphere (section 3 5) and reduce the present discrepancy. 

The vertical and latitudinal magnitude of the cooling, and likewise the location of the 
warming region above the cooling in the lower stratosphere, are very sensitive to the vertical 
profile of ozone depletion imposed in the model. The models invariably locate the cooling at 
exactly the region of the imposed ozone loss, with a warming immediately above it at the higher 
latitudes. Thus, any shift of the altitude extent of ozone depletion in the model has the potential 
to shift the peak cooling, and thus alter the vertical profile of the computed cooling trend. In 
turn, this affects the quantitative inferences about the consistency between computed and ob- 
served temperature trends in the middle and upper stratosphere. 

4.3 Effects in upper troposphere 

The cooling of the lower stratosphere in response to ozone losses there also leads to a cool- 
ing below the ozone loss region in the global upper troposphere (-10-14 km.), in part due to 
reduced infrared emission from the stratosphere (Ramanathan and Dickinson, 1979). The cool- 
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ing of the lower stratosphere also impacts the forcing of the troposphere (WMO: 1992, 1995; 
Chapter 10). Specifically, an ozone-induced cooling of the lower stratosphere implies a reduc- 
tion in the longwave radiative emission from the stratosphere into the troposphere. This leads 
to a negative radiative forcing of the surface-troposphere system (Hansen et al., 1997a). In gen- 
eral, GCMs simulating the effects of stratospheric ozone loss indicate a cooling of the upper 
troposphere (see Figures 25, 27, 28). 

All modeling efforts that have considered appropriate changes in trace gases yield a cooling 
of the tropopause region (e.g., Hansen et al., 1995, see also Figure 27; Vinnikov et al., 1996; 
Santer et al., 1996; Hansen et al., 1997a). Here, again, the lower stratospheric ozone depletion 
effects dominate over all other gases including C0 2 . Further, the vertical profile of changes in 
ozone critically determines the profile of the temperature change in the upper troposphere. 

A difficulty arises in comparing simulated upper tropospheric temperature changes due to 
ozone losses with observations over the past decade as the ‘signal’ from models (see Figures 
27, 28) is quite small - less so than that in the lower stratosphere (say, at 50 hPa or ~20 km.). It 
may be noted that the modeled upper tropospheric changes, even for as a large a forcing as a 
doubling of C0 2 , are generally small compared to those in the lower troposphere or higher in 
the stratosphere (Vinnikov et al., 1996). It is difficult at present to interpret the small change in 
the upper troposphere temperatures inferred from observations in terms of a change originating 
primarily from the stratosphere (e g., lower stratospheric ozone loss) or one arising in the tro- 
posphere (e.g., tropospheric ozone chemistry). 

In general, the temperature changes near the tropopause have the potential for affecting the 
water vapor distribution and cloudiness, especially cirrus clouds, which in turn can cause feed- 
back effects in the climate system. One example available from observations of the effects in 
the upper troposphere due to a stratospheric perturbation concerns the aftermath of the Pinatu- 
bo eruption. Lidar and other observations (Sassen et al., 1995) suggest that the aerosols from 
this volcanic eruption affected cirrus (ice) cloud formations near the tropopause. 

The interactions between trace species’ changes and the associated chemistry, microphys- 
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ics, radiation and dynamics remain to be more thoroughly explored, both from modeling and 
observational perspectives, taking into account the three-dimensional nature of the problem 
(Holton et al., 1995). Meteorological parameters such as tropopause location (Schubert and 
Munteanu, 1988; Hoinka et al., 1998a), temperature and geopotential height (Ohring and 
Muench, 1960; Spankuch and Schulz, 1995) tend to be highly correlated to total ozone. They 
can, therefore, be expected to exhibit changes owing to trends in ozone. At this time, such 
changes are not well-documented (section 2.7). There is growing evidence of a connection be- 
tween stratospheric variability and trends with those in the troposphere (e.g., Kodera et al., 
1994; Kitoh et al., 1996; Perlwitz and Graf, 1995). As noted earlier, Thompson and Wallace 
(1998) show a connection between the trend in winter-mean NH surface temperature and pres- 
sure and the trend in polar vortex strength. Such linkages inferred from observations emphasize 
that the radiative-chemical-dynamical interactions, which couple the stratospheric and tropo- 
spheric processes, need to be fully accounted for in order to attribute the observed temperature 
variations and trends to well-defined mechanisms that include changes in trace species. 
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5. SUMMARY 


I. Observations 

(a) Datasets available for analyzing stratospheric temperature trends comprise measurements by 
radiosonde ( 1 940s-present), satellite ( 1 979-present), lidar ( 1 979-present) and rocketsonde (periods 
varying with location, but most terminating by -mid- 1990s), meteorological analyses based on ra- 
diosonde and/or satellite data, and products based on assimilating observations using a general cir- 
culation model. 

(b) The temporary global annual-mean lower stratospheric (-50-100 hPa) warming (peak value 

- 1 K) associated with the aerosols from the Pinatubo volcanic eruption (see WMO, 1992 and 1995), 
which lasted up to about 1993, has now given way to a relatively colder stratosphere. 

(c) Radiosonde and satellite data indicate a cooling trend of the global, annual-mean lower strato- 
sphere since - 1980. Over the period 1979-1994, the trend is -0.6K/ decade. For the period prior 
to 1980, the radiosonde data exhibits a substantially weaker long-term cooling trend. 

(d) Over the period 1979-1994 there is an annual-mean cooling of the Northern Hemisphere mid- 
latitude lower stratosphere (-0.75K/ decade at 30-60N). This trend is coherent amongst the various 
datasets with regard to the magnitude and statistical significance. Over the longer period 1966- 
1994, the available datasets indicate an annual-mean cooling at 30-60N of -0.3K/ decade. 

(e) In the -15-45 degree latitude belt of the Southern Hemisphere, the radiosonde record indicates 
an annual-mean cooling of the lower stratosphere up to -0.5-1K/ decade over the period 1979- 
1994. The satellite record also indicates a cooling of the lower stratosphere in this latitude belt, but 
this is statistically significant only between about November and April. 

(0 Substantial cooling (-3-4K/ decade) is observed in the polar lower stratosphere during late win- 
ter/spri ngti me in both hemispheres. An approximate decadal-scale cooling trend is evident in the 
Antarctic since about the early 1 980’s, and in the Arctic since about the early 1990’s. However, the 
dynamical variability is large in these regions, particularly in the Arctic, and this introduces diffi- 
culties in establishing a high statistical significance of the trends. 

(g) A cooling of the upper stratosphere (p<3 hPa; Z>40 km.) is apparent over the 60N-60S region 
from the annual-mean SSU satellite data over the 1979-1994 period (up to -3K/ decade near 50 
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km). There is a slight minimum in cooling in the middle stratosphere (-30-40 km.) between the 
maxima in the lower and upper stratosphere. 

(h) Lidar and rocket data available from specific sites generally show a cooling over most of the 
middle and upper stratosphere (-30-50 km.) of 1 to 2K / decade since -1970, with the magnitude 
increasing with altitude. The influence of the 1 1-year solar cycle is relatively large (>1K) at these 
altitudes (>30 km.) 

(i) The vertical profile of the annual-mean stratospheric temperature change observed in the North- 
ern Hemisphere midlatitude (45N) over the 1979-1994 period is robust among the different 
datasets. Figure 32 illustrates the overall mean vertical profile of the trend and uncertainty at 45N, 
taking into account all of the datasets and accounting for the uncertainties of the individual mea- 
surements. The overall trend consists of a -0.75K/ decade cooling of the -20-35 km. region, with 
the cooling trend increasing with height above (-2.5 YJ decade at 50 km.). The vertical profile of 
cooling, and especially the large upper stratospheric cooling is consistent with the global plots in 
WMO (1988 (e.g., Figure 6.17); 1990 (e.g., Figure 2.4-5)) constructed from shorter data records. 


EL Model Results and Model-Observation Comparisons 

(a) Model simulations based on the known changes in the stratospheric concentrations of various 
radiatively-active species indicate that the depletion of lower stratospheric ozone is the dominant 
factor in the explanation of the observed global-mean lower stratospheric cooling trend (-0.5-0.6 
YJ decade) for the penod 1 979- 1 990. The contribution to this trend from increases in well-mixed 
greenhouse gases is estimated to be less than one-fourth that due to ozone loss. 

(b) Model simulations indicate that ozone depletion is an important causal factor in the latitude- 
month pattern of the decadal (1979-1990) lower stratosphere cooling. The simulated lower strato- 
sphere in Northern and Southern Hemisphere midlatitudes, and in the Antarctic springtime, gener- 
ally exhibits a statistically significant cooling trend over this period consistent with observations. 

(c) The Fixed Dynamical Heating or FDH (equivalently, the pure radiative response) calculations 
yield a mid-to-high latitude annual-mean cooling that is approximately consistent with a General 
Circulation Model’s (GCM) radiative-dynamical response (Figure 28); however, changes in circu- 
lation simulated by the GCM cause an additional cooling in the tropics, besides affecting the me- 
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ridional pattern of the temperature decrease. 

(d) FDH model results indicate that both well-mixed greenhouse gases and ozone changes are im- 
portant contributors to the cooling in the middle and upper stratosphere; however, the computed 
upper stratospheric cooling is smaller than the observed decadal trend. Increased water vapor in the 
lower-to-upper stratosphere domain could also be an important contributor to the cooling; however, 
decadal-scale global stratospheric water vapor trends have not yet been determined. 

(e) Model simulations of the response to the observed global lower stratospheric ozone loss in mid- 
to-high latitudes suggest a radiative-dynamical feedback leading to a warming of the middle and 
upper stratospheric regions, especially during springtime; however, while the modeled warming is 
large and can be statistically significant during the Antarctic spring, it is not statistically significant 
during the Arctic spring. Antarctic radiosonde observations indicate a statistically significant 
warming trend in Spring at -30 hPa (24 km.) and extending possibly to even higher altitudes; this 
region lies above a domain of strong cooling that is approximately co-located with the altitude of 
observed ozone depletion. 

(0 There is little evidence to suggest that tropospheric climate changes (e.g., induced by green- 
house gas increases in troposphere) and sea-surface temperature variations have been dominant 
factors in the global-mean stratospheric temperature trend over the 1979-1994 period. The effect 
of potential shifts in atmospheric circulation patterns upon the decadal trends in global stratospher- 
ic temperatures remains to be determined. 
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TABLE 1 


Zonal temperature time series made available to and considered by SPARC-STTA. Angell, Oort, 
Russia, Raob and Berlin are different radiosonde datasets. MSU and Nash are satellite instruments 
while lidar data is from OHP (France). CPC, Reanal, UKMO/ SSUANAL and GSFC are analyzed 
datasets. For the MSU and Nash satellite data, the approximate peak levels ‘sensed’ are listed. Ref- 
erences to earlier versions of the datasets are also listed. See section 2.1 for details. 


Dataset 


Period 


Location Monthly 


. Levels (hPa) 


Angell 

[Angell, 1988] 


1958-1994 8 bands 

4 bands 


3-monthly 100-50 

3-monthly 50, 30, 20, 10 


Oort 1958-1989 85S-85N 

[Oort and Liu, 1993] 


monthly 100, 50 


Russia 


1959-1994 

1961-1994 


[Koshelkov and Zakharov, 1998] 


7 ON, 80N 
70N, 80N 


monthly 100 

monthly 50 


UK Raob (or Raob) 1961-1994 87.5S-87 5N monthly 100,50,30,20 

[Parker and Cox, 1995] 


Berlin 1965-1994 10-90N monthly 100,50,30 

[Labitzke and van Loon, 1994] 


Lidar 1979-1994 44N, 6E monthly 10,5,2,1,0.4 

[Hauchecome et al., 1991] 


MSU 1979-1994 85S-85N 

[Spencer and Christy, 1993] 


Nash 1979-1994 75S-75N 

[Nash and Forrester, 1986] 


monthly 90 


monthly 50, 20, 15, 6, 5 

2, 1.5, 0.5 
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CPC 

1979-1994 

85S-85N 

monthly 

[Gelman et al., 1994] 

1964-1978 

20N-85N 

monthly 

Reanal 

[Kalnay et al., 1996] 

1979-1994 

85S-85N 

monthly 

GSFC 

[Schubert et al., 1993] 

1979-1994 

90S-90N 

monthly 

UKMO/ SSUANAL 
[Bailey et al., 1993] 

1979-1994 

90S-90N 

monthly 
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70, 50, 30, 10, 
5,2, 1 
50, 30, 10 


100, 70, 50, 
30, 10 


100, 70, 50, 
30, 20 


50, 20, 10, 5, 

2 , 1 



TABLE 2 


Rocketsonde locations and periods of coverage utilized for the present assessment [based on 
Dunkerton et al., 1998; Golitsyn et ah, 1996; Keckhut et ah, 1998; Kokin and Lysenko, 1994; Ly- 
senko et ah, 1997; and Komuro, 1989 (updated)]. 


Station 

Latitude: Longitude 

Period 

Heiss Island 

81N; 58E 

1964-1994 

Volgograd 

49N; 44E 

1965-1994 

Balkhash 

47N; 75E 

1973-1992 

Ryon 

39N, 141.SE 

1970-present 

Wallops Island 

37.5N; 76W 

1965-1990 

Point Mugu 

34N; 1 19W 

1965-1991 

Cape Kennedy 

28N; 80W 

1965-1993 

Barking Sands 

22N; 160W 

1969-1991 

Antigua 

17N; 61W 

1969-1991 

Thumba 

08N; 77E 

1971-1993 

Kwajalein 

09N; 167E 

1969-1990 

Ascension Island 

08S; 14W 

1965-1993 

Molodezhnaya 

68S; 46E 

1969-1994 
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TABLE 3 


Lower stratospheric temperature trends from published studies of radiosonde data. 


Reference 

Data 

Period 

Level or 
Layer 

Region 

Trend (K/ 
decade) 

Comments 

Angell, 1991a 

1972-198 

9 

50hPa,20km 

8S - 55N 

-0.5±0 3 

12 radiosonde stations. 
Data were adjusted for 
the El Chichon influence. 
Greatest cooling in 
winter. 

30 hPa,24km 


20hPa,27km 

mm 

10hPa,31 km 

-0.3±0 3 

Angell, 1991b 

1970-198 

8 

100-50 hPa 
layer 

NH 

-0.2 

Trends reported here are 
based on Angell’s 
presentation of 
temperature differences 
between the periods 
1980-88 and 1970- 78. 

SH 

-0.5 

Koshelkov and 

Zakharov, 

1998 

1965-199 

4 

50 and 100 
hPa levels 

65N - 83N 

0 to -1 (-0.5?) 

Trends estimated from 
graphed results. 
Significant (95 % 
confidence level) trends 
only from May to 
October 

1979-199 

4 

0 to -1 

Trends estimated from 
graphed results. 
Significant (95 % 
confidence level) trends 
only from June to August 

Labitzke and 
van Loon, 1995 

1965-199 

3 

100 hPa 

NH (10-90N) 

0 to -0.5 

Significant (95% 
confidence level) trends 
between 60 and 80N 

50 hPa 

NH (10-90N) 

0 to -0.5 

Significant (95 % 
confidence level) trends 
between 20 and 90N 

30 hPa 

NH (10-90N) 

0 to -0.5 

Significant (95% 
confidence level) trends 
between 30 and 60N 

1979-199 

3 

100 hPa 

NH (10-90N) 

-0.2 to -0.8 

Significant (95 % 
confidence level) trends 
only at about 40N 

50 hPa 

NH (10-90N) 

0 to -0.9 

Significant (95 % 
confidence level) trends 
between 35 and 50N 

30 hPa 

NH (10-90N) 

0 to -1.0 

Significant (95 % 
confidence level) trends 
between 30 and 50N 
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Miller et al., 
1992 

1964-198 

6 

150-30 hPa 
levels 

Global (62 
stations, as 
Angell) 

-0.2 to -0.4, 
depending on 
level 

McCormack 
and Hood, 1994 

1979-199 

0 

100 hPa 

NH 

0 to -4.5 

Oort and Liu, 
1993 

1963-198 

8 

100-50 hPa 
layer 

NH 

E 

.38+0.14 



SH 

-0.43+0 16 


Globe 

-0.40+0 12 

1959-198 

8 

NH 

1 

.40+0 10 

Pawson and 
Naujokat, 1997 

1965-199 

6 

50 hPa 

NH (10-90N) 

-1.90(minimum 

) 

+1.67 

(maximum) 



30 hPa 


-1.85 

(minimum) 

+4.07 

(maximum) 

Parker et al., 
1997 

1965-199 

6 

150-30 hPa 
layer 

NH 

-0.27 



SH 

-0.44 

1979-199 

6 


NH 

-0.63 


SH 

-0.73 

Reid et al., 1989 

1966-198 

2 

100-15 hPa 
levels 

tropics 

-1.2 to +0.8 

Taalas and 
Kyro, 1992 

1965-198 

8 

50 hPa 

Sodankyla, 

Finland 

station 

-0.2±1 6 
(annual) 
-1.6±0 08 
(January) 
+1.5±0 05 
(Apnl) 

Taalas and 
Kyro, 1994 

1965-199 

2 

50 hPa 

Sodankyla, 

Finland 

station 

-0.16 


Four stations’ data were 
adjusted for level shifts. 


Data from the FUB 
analyses. Significant 
trends only in winter at 
30-45N. 



Decrease in daily 
minimum temperature, 
increase in daily 
maximum, for winter 
season. Increase in the 
area of T<195K and 
T<192K, although the 
early data for T<192 are 
questionable. 


Radiosonde data 
weighted to correspond 
with MSU4. 

Adjustments made to 
Australasia data for 
1979- 1996 . Trends 
significant at the 95% 
confidence level or better. 


Five radiosonde stations. 
Trends vary by station 
and level. 


100, 70 and 30 hPa 
showed similar results 


Also, an increase in the 
number of observations 
of T<195K in winter. 
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TABLE 4 


Latitude bands where the observed 50hPa, 100 hPa and satellite temperature trends (1979-1994) 
from various data sources (see Table 1 and Figure 5) are statistically significant at the 2-sigma level. 
SH and NH denote southern and northern hemisphere, respectively. No entry denotes either no data 
or no statistically significant latitude belt in that hemisphere. 

50 hPa 


Data 


Latitude Band 


SH 


NH 


Berlin 

- 

30-55N 

GSFC 

- 

42-58N 

Reanal 

37.5-25S 

32.5-55N 

Raob 

12.5S 

32.5-62. 5N 

CPC 

30S 

10N 

Angell 

- 

50N 

Russia 

100 hPa 


Data 

Latitude 

; Band 


SH 

NH 

Berlin 

_ 

0-5N; 35-80N 

GSFC 

- 

38-60N 

Reanal 

- 

37. 5-57. 5N 

Raob 

47.5S 

12.5N; 37.5-72.5N 

Russia 

Satellite 

7 ON 

Data 

Latitude Band 


SH 


NH 


Nash 

MSU 


30-60N 
28. 75-63. 75N 
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TABLE 5 


Latitude bands where the 1 966- 1 994 observed 50hPa temperature trends (1979-1 994) from various 
data sources (see Table 1 and Figure 6) are statistically significant at the 2-sigma level. SH and NH 
denote southern and northern hemisphere, respectively. No entry denotes either no data or no sta- 
tistically significant latitude belt in that hemisphere. 


50 hPa 


Data Latitude Band 



SH 

NH 

Berlin 

- 

20-70N 

CPC 

- 

10-20N; 50-60N 

Angell 

- 

10, 20, 50N 

Raob 

37.5- 32.5S; 

22.5- 17.5S 

17.5N; 32.5N; 42.5-72.5N 

Oort 

84-56S; 30-28S; 
24-22S; 18-16S 

20-30N; 34N; 40N; 46N; 
52-68N 

Russia 

_ 

70N 
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FIGURE CAPTIONS 


Figure 1 . Altitude range of the signals ‘sensed’ by the various thermal infrared channels of the 
Stratospheric Sounding Unit instrument (SSU), and that ‘sensed’ by channel 4 of the Microwave 
Sounding Unit (MSU). (Left) Nadir and (right) off-nadir channel weighting functions. 

Figure 2. Time series of global temperature anomalies from the overlap-adjusted SSU data. These 
data measure the thermal structure over thick layers (-10-15 km.) of the stratosphere; the label 
in each panel indicates the approximate altitude of the weighting function maximum. 

Figure 3. (Top) Global and hemispheric averages of annual anomalies of 100-50 hPa layer-mean 
virtual temperature from 63 radiosonde stations (Angell, 1988, updated; Halpert and Bell, 
1997). (Bottom) Same as (top), except from MSU channel 4 (see section 1 for stratospheric 
altitudes ‘sensed’); the thick, solid line denotes global-mean, the thin, solid line NH mean and 
the dashed line the SH mean. [Updated from Randel and Cobb, 1994]. 

Figure 4. 1 955-1 997 time series of July 30 hPa temperatures at high northern latitude (80N). Trends 
over the 1955-1977, 1979-1997 and 1955- 1997 time periods are -0.0 1,-1. 41 and -0.5 K/ decade, 
respectively. Of these, only the 1979-1997 trend is statistically significant. [Updated from 
Labitzke and van Loon, 1995]. 

Figure 5. Zonal-mean decadal temperature trends for the 1979-1994 period, as obtained from 
different datasets. These consist of radiosonde (Angell, Berlin, UKMO/ Raob and Russia/ Raob) 
and satellite observations (MSU and Nash), and analyses datasets (CPC, GSFC and Reanal). 
See Table 1 , and sections 2 1 and 2.3 2 for details, (a) denotes 50 hPa trends, (b) denotes 100 
hPa trends, and (c) denotes trends observed by the satellites for the altitude range ‘sensed’ which 
includes the lower stratosphere (see Figure 1). Latitude bands where the trends are statistically 
significant at the 2-sigma level are listed in Table 4. [Data courtesy of S PARC-Stratospheric 
Temperature Trends Assessment project]. 

Figure 6. Zonal-mean decadal temperature trends at 50 hPa over the 1966-1994 period from 
different datasets (see Table 1 , and sections 2. 1 and 2.3.2 for details). Latitude bands where the 
trends are statistically significant at the 2-sigma level are listed in Table 5. [Data courtesy of 
SPARC-Stratospheric Temperature Trends Assessment project]. 

Figure 7. Zonal, annual-mean decadal temperature trends versus altitude for the 1979-1994 period, 
as obtained from the Nash satellite (including MSU and SSU channels) retrievals (top) with 
volcanic periods included and (bottom) with volcanic periods omitted. Shaded area denotes 
significance at the 2-sigma level. [Data courtesy of SPARC-Stratospheric Temperature Trends 
Assessment project]. 

Figure 8. Trend (K/ year) from the all-seasons lidar record at Haute Provence, France (44N, 6E) 
over the period 1979-1998. The 2-sigma uncertainties are also indicated. [Updated from 
Keckhut et al., 1995]. 

Figure 9 (top, middle and bottom panels). Vertical profile of the zonal, annual-mean decadal 
stratospheric temperature trend over the 1 979- 1 994 period at 45N from different datasets (Table 
1 ) for the 1979-1 994 period. Horizontal bars denote statistical significance at the 2-sigma level 
while vertical bars denote the approximate altitude range ‘sensed’ by the MSU, and by the 
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different SSU satellite channels (see Fig. 1). [Data courtesy of SPARC-Stratospheric 
Temperature Trends Assessment project] 

Figure 10. (A) Trends (K/ year) based on the Former Soviet Union (FSU) rockets (25-70 km.) since 
the mid-1960s (see Table 5-2) [Adapted from Golitsyn et al., 1996]. (B) Trends (K / year) based 
on the same FSU rocket dataset. [Adapted from Lysenko et al., 1997], 

Figure 1 1 . (A) Trend (K/ year) compiled from US rocket sites in the 8S-34N belt (see Table 2) over 
the 1969-1993 (Keckhut et al., 1998). (B) Trend (K/ year) from the Japanese rocket station at 
Ryori over the period 1970-present. [Updated from Komuro, 1989]. 

Figure 12 (top, middle and bottom panels). Vertical profile of the zonal, annual-mean decadal 
stratospheric temperature trend over the 1979-1994 period at 28-38 N from different datasets 
(Tables 1 and 2) for the 1 979- 1 994 period. Horizontal bars denote statistical significance at the 
2-sigma level while vertical bars denote the approximate altitude range ‘sensed’ by the MSU 
and by the different SSU satellite channels (see Figure 1). [Data courtesy of SPARC-Strato- 
spheric Temperature Trends Assessment project]. 

Figure 13. Latitude-time section of trends in MSU Channel 4 temperature over January 1979-May 
1998. Dashed lines indicate cooling trends. Contour interval is 0.5 K/decade, and shading in- 
dicates the trends are significant at the 5% level. [Updated from Randel and Cobb, 1994]. 

Figure 14. (Top) Time series of 100 hPa zonal-mean temperatures at 80N in March from NCEP 
reanalyses (solid lines), together with data from MSU. The smooth curve through the reanalysis 
data indicates the decadal-scale variation. (Bottom) Time series of 100 hPa temperature in No- 
vember at Halley Bay, Antarctica, from radiosonde data (solid lines) and NCEP reanalyses 
(dashed lines). The smooth curve through radiosnde data indicates the decadal-scale variation. 
[Randel and Wu, 1998]. 

Figure 15. Altitude-time section of temperature differences over Antarctica, 1987-96 minus 1970- 
79, derived from an average of seven radiosonde stations. Dashed lines indicate cooling trends. 
Contour interval is IK. Shading denotes that the temperature differences are significantly dif- 
ferent from natural variability (Randel and Wu, 1998). 

Figure 16. Seasonal rocketsonde trends in K/ year, using FSU rocket data from 1972 to 1990 at 5 
sites (see Table 2); winter (solid line), spring (dot-dash line), summer (dotted line) and autumn 
(dashed line). [Adapted from Kokin and Lysenko, 1994]. 

Figure 17. Monthly-mean 50hPa raw temperatures measured by radiosondes at Tahiti. Several 
instrument changes occurred during this data period, but the 1976 change from Mesural FMO 
1943B, with a bimetal temperature sensor, to Mesural FMO 1944C, with a thermistor, had the 
most obvious effect on the time series. [Gaffen, 1996; see also Gaffen, 1994]. 

Figure 18. Time series of deseasonahzed temperature anomalies at the equator for NCEP data at 1 
hPa (top panel) and SSU 27 data (bottom panel). Vertical dashed lines in top panel indicate 
changes in the operational satellites. Refer to Figures 1 and 2 concerning SSU 27. 
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Figure 19. Latitude-altitude structure of the solar signal in the Nash satellite dataset, derived from 
regression analyses. Shaded areas show 95% confidence intervals for the regression estimates. 
[Data courtesy of SPARC-Stratospheric Temperature Trends Assessment project]. 

Figure 20. Mean temperature trend for 6 US rocket stations indicating the influence of the solar 
cycle (Dunkerton et al., 1998). 

Figure 21. Change in tropopause height compared with the change in the 5-km. temperature at Ho- 
henpeissenberg. The 5-km. temperature (dark line) shows a trend of 0.7K7 decade while the 
tropopause height (grey line) shows a trend of 150m/ decade (Steinbrecht et al., 1998). 

Figure 22. Computed global- and annual-mean temperature changes in the -50-100 hPa (-16-21 
km.) lower stratospheric region. The result for the ozone depletion corresponds to the 1979 to 
1990 losses, and is computed using a GCM. The results for C0 2 alone, and for all the well- 
mixed greenhouse gases together (viz., C0 2 , CH 4 , N 2 0 and CFCs) are computed using a ID 
RCM for three different periods. [Reprinted by permission from Nature (Ramaswamy et al., 
Nature, 382, 616-618, 1996) Copyright (1996) Macmillan Magazines Ltd.]. 

Figure 23. Vertical profile (-15-35 km.) of the FDH-computed temperature change during the 
1980s at 40N due to changes in the concentrations of various stratospheric species. “Green- 
house gases” denotes well-mixed greenhouse gases; the “SBUV” ozone refers to calculation 
which considers the column loss using SBUV satellite trends; the “SAGE-B” calculation uses 
the vertical profile of ozone change observed by the SAGE satellite; “water vapor” refers to the 
decadal change in water vapor reported by Oltmans and Hoffman (1995). [Adapted from For- 
ster and Shine, 1997]. 

Figure 24. Vertical profile of the temperature change from the lower to upper stratosphere resulting 
due to the 1979-1990 changes in the well-mixed greenhouse gases. [The results are from the 
FDH model used in Ramaswamy et al., 1992]. 

Figure 25. GCM-computed temperature change in the atmosphere due to a total loss of ozone be- 
tween 70 and 250 hPa. Results are shown without (dashed line) and with (solid line) consider- 
ation of cloud feedbacks. [Adapted from Hansen et al., 1997]. 

Figure 26. Modeled evolution of temperature changes in the lower stratosphere (80S) using two 
different assumptions. FDH denotes the usual Fixed Dynamical Heating while SEFDH denotes 
a Seasonally Evolving Fixed Dynamical Heating. [Adapted from Forster et al., 1997]. 

Figure 27. Comparison of GCM-simulated vertical profile of temperature change (K) over the 
1979-1990 period with radiosonde and MSU satellite observations. The simulations are per- 
formed considering well-mixed greenhouse gas increases, and with or without consideration of 
the lower stratospheric ozone depletion. [Adapted from Hansen et al., 1995]. 

Figure 28. Annual-mean stratospheric ozone loss profile (top panel), and the corresponding tem- 
perature changes, as obtained using a FDH model (middle panel) and GCM (bottom panel). 
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[Adapted from the model simulations of Ramaswamy et al. (1992 and 1996)]. 

Figure 29. Zonally-averaged, monthly mean, lower stratospheric temperature change 1979-1991: 
(a) as simulated by the general circulation model (90S to 90N) due to the observed global ozone 
depletion, and (b) as inferred (Randel and Cobb, 1994) from satellite observations (82.5S to 
82. 5N). Shaded areas show statistical significance at the 95% confidence level. [Reprinted by 
permission from Nature (Ramaswamy et al., Nature, 382, 616-618, 1996) Copyright (1996) 
Macmillan Magazines Ltd.]. 

Figure 30. FDH-computed global stratospheric temperature change during the 1980s versus lati- 
tude due to (a) SBUV column depletion with ozone loss between local tropopause and 7 km. 
above it, and (b) assumption of the SAGE vertical profile. [Adapted from Forster and Shine, 
1997]. 

Figure 3 1 . Comparison of low latitude (5N) temperature response due to solar cycle variation, as 
computed by two models (Brasseur, 1993, and McCormack and Hood, 1996), with SSU satel- 
lite observations (thick, dashed line), and rocket data (Dunkerton et al., 1998) for the altitude 
range -28-55 km. (vertical bar). McCormack and Hood perform two FDH calculations using 
the range of the observed ozone values; their model results (solid curves) are quite similar to 
that of Brasseur (1993, thin dashed curve). [Adapted from McCormack and Hood, 1996]. 

Figure 32. Summary Figure illustrating the overall mean vertical profile of temperature trend over 
the 1979-1994 period in the stratosphere at45N, as compiled using radiosonde, satellite and 
analyzed datasets (section 2 3.3) The combined uncertainty due to the measurements at the 
various altitudes is also indicated. [Data courtesy of SPARC-Stratosphenc Temperature Trends 
Assessment project]. 
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